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ABSTRACT 


The PHALANX Close-In-Weapons System provides the last layer of defense 
against an incoming anti-ship missile. Current dis|>ersion data shows that the 
rotating six-barrel system firing at 75 shots/second gives variable dispersion 
greater than 2 mradian. This results in a 4 m spread at its maximum 2000 m 
range. Computer simulations have shown that a fast missile can penetrate to 
within a couple hundred yards of the ship, where even if it was destroyed, the 
resulting ballistic debris would hit the ship in a shotgun blast pattern. 


In order to investigate the physics behind the dispersion and test future 
design changes, a detailed finite element model of the PHALANX gun was 
developed. This is the first finite element model of the gun ever. The normal 
modes of a single barrel have been matched favorably with those of a real barrel. 
This model has been used to calculate the normal modes of vibration of the whole 
gun in order to predict the motion of the complex system. Barrel tip displacement, 
which plays a critical role in the dispersion pattern, was determined after applying 
the periodic firing excitation force and damping. Its motion corresponds well to the 
2 mradian dispersion. In addition, its first mode of vibration at 14 Hz appears to 
be significant in the analysis of the dispersion data. This analysis has focused 
attention on certain components of the gun, such as the double angular contact 
bearing, which may be a leading cause of the dispersion pattern. Finally, this gun 
model can also support design changes such as the addition of barrel restraints or 


the proposed new barrels and their effect on dispersion predicted prior to any 


actual hardware tests or procurement orders. 
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I. INTRODUCTION 


The PHALANX Close-In-Weapon-System (CIWS) is considered to be a ship’s 
last layer of defense against incoming missiles. The PHALANX CIWS is shown in 
Figure 1-1. Nearly every ship of the current 465 ship U. S. fleet has at least one of 
these gun mounts. For some ships, it is the only means of missile defense. 
Although it has been in the fleet for over a decade, it has never been tested in 
combat, so its effectiveness has largely been measured by either missile range tests 
or computer simiilations. Within the past flve years, there has been documented 
concern that this defensive weapon system may not be that effective, especially in 
defeating sophisticated missiles that are currently available to Third World 
countries. The problem of its effectiveness is attributed almost entirely to that of 
the gun itself and not its Are control system [Ref. 1]. At the maximum range of the 
PHALANX gun, the incoming missile or patrd boat presents a small cross-section. 
The angular dispersion of bullets leaving the gun can dimmish the gun’s ability to 
hit targets, especially as the target’s size decreases, becomes more maneuverable, 
or increases its velocity. 

As part of a class project, supervised by Professor W. B. Colson, the problems 
associated with the gxm were illustrated in some of the best case scenarios. 
Specifically, a computer simulation was used to study a PHALANX engagement of 
an incoming anti-ship missile. The results of this class project, discussed in 
Chapter II, provided the motivation for this thesis. 

The purpose of this study has been to use a sophisticated computer code to 
create a finite element model of the PHALANX gun in order to better imderstand 
its complex dyneunic behavior. The Naval Postgraduate School has access to one of 
the most advanced finite element codes available in private industry today, called 
I-DEAS (Integrated Design Engineering Analysis) which is developed and 


1 




Figure 1*1. The PHALANX Cloee-ln Weapon System [Ref. 2], 


supported by SDRC (Structural Dynamics Research Corporation). The finite 
element model was used to predict barrel tip displacement through such causes as 
barrel oscillations and the coupling of the barrels to the gun mount through its 
bearings. In addition, the finite element model provides an important tool for the 
evaluation of future design changes of the PHALANX gun system. Chapter III 
discusses the benefits of using this type of vomputer code and the theory behind its 
development. 

Chapter IV then discvisses how raw data from blueprints and technical 
discussions with engineers familiar with the PHALANX CIWS was molded into a 
model of the gxm. Field trips to Hughes Missile Systems Company, Pomona, 







California, and the USS ARKANSAS provided further modeling information. 
Valuable finite element modeling experience by Mike Hatch, an expert on rotating 
computer disk systems >ind Professor J. H. Gordis, a former SDRC mechanical 
vibration theorist, ir^atly enhanced the effectiveness of the model. 

In Chapter V, the normal modes of vibration are calculated for the PHALANX 
gun. Solving for the normal modes is an important and insig^itful step into 
understanding the dynamic response of any structure. Determining which normal 
modes of a structure are excited after applying particxilar excitation forces and 
setting boundary conditions ultimately determines how that structure will behave. 
In analogy, any color of lig^t can be made up of the three primary colors, red, blue, 
and yellow. The response of any system can be defined in terms of its normal 
modes. The dynamic response package of I-DEAS uses modal degrees of freedom 
instead of physical degrees of fineedom, in order to greatly reduce the dynamic 
response solving time. 

Chapter VI discusses how these normal modes are used in a dynamic 
response analysis. Damping is defined for the gtm system and the primary driving 
force of the gun while firing is applied. The resulting motion of the barrel tip, 
which gives a good indication of the dispersion pattern, is then displayed for each 
of the three translational directions. More importantly, the firequency response 
functions are also calculated and are used to reveal which normal modes are being 
excited and to what extent. Such an analysis may ultimately lead to tuning a 
problem mode so that its excitation is limited, and the resulting barrel tip motion 
is correspondingly reduced. 

Recommendations for follow-on modifications of the model are listed in 
Chapter VII. Although the current model has been correlated in a number of ways 
with the real gun, it is important to refine the model further so as to enhance its 
accuracy while reducing its total number of degrees of fireedom. By reducing the 
model’s degrees of freedom, the computer solve time is shortened which can lead to 
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a smoother animation picture, thus providing better insight into its dynamic 
behavior. Also, if the proposed muzzle restraint and isogrid is added, it will be 
important to keep the model’s complexity within a manageable budget. 

Although the primary purpose of this study was to build and demonstrate the 
practical value of a finite element model of the gun, some preliminary analysis of it 
was performed and the resiilts are discussed in Chapter VIII. Of key interest is 
the first flexible mode of vibration at 14 Hz and the double angular contact bearing 
and the effects they may have on the dispersion pattern. 

The purpose of this study has been multi-faceted. While this computer model 
will be important in predicting the observed dispersion pattern, it can also serve as 
a program manager’s tool in helping make future design decisions. Components 
can be swapped out with new ones, and their effect on the dynamic behavior of the 
gun can be predicted long before any hardware tests are actually performed. In 
summary, a tool of the 1990’s is being used to improve 1950’8 technology. 
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II. MOTIVATION FOR STUDYING THE PHALANX GUN 


A. DESCRIPTION OF THE PHALANX SYSTEM 

The PHALANX Close-In Weapon System (CIWS) is a ship-based, stand alone 
gun mount designed to disable an incoming missile or, in a recent application, 
engaging fast patrol boats. It has been actively deployed for over a decade in 
virtually every ship of the U. S. fleet, and has had to undergo continual review and 
various modifications to match the ever ’’fluid" world threat. In that past decade, 
the missile it was initially designed to destroy has become faster, smaller, smarter 
and more maneuverable. Furthermore, the U. S. fleet has also made a transition 
from operating in the open, "blue water," sea to the coastal, "brown water," 
environment. With this change in operating areas, small speed boats, which were 
considered an insignificant challenge in blue water, now have exposed a weakness 
in surface combatants. Rather than building a new weapon system to deal with 
this threat, the U. S. Navy has decided to adapt the PHALANX CIWS into an 
anti-patrol boat mode. These two growing threats, with their decreasing radar 
cross-section, and increasing ability to pop-up with little notice, requires the 
PHALANX to move beyond its original design envelope. The primary area of 
concern with the PHALANX gun is to tighten its dispersion pattern to put more 
rounds on target, sooner and hence destroy the target at a greater range from the 
ship. This is the focus of this study. 

Before discussing the natiure of the problem, it is important to have an overall 
understanding of the gun system. Figure 2-1 shows a component breakdown of the 
PHALANX weapon system. Starting at the top, the dome structure contains an air 
search radar which provides the initial threat detection. It has a maximum 
effective range for detection greater than 4500 m (5000 yds). It provides 360 
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degrees, horizon to azimuth search coverage. Once a target has been identified by 
the search radar and the computer determines that it meets the criteria of a threat 
to the ship, the gxm moimt slews to the target bearing and transfers tracking of it 
to the fire control radar. This radar is mounted below the search radar, on the 
side of the cylindrical structure, and is physically boresighted to the gun boresight. 
It is a high resolution radar whose purpose is to track the target as well as the 
outgoing rounds. The fire control computer sends the necessary servo information 
to keep this tracking radar, and hence the gun, locked on target [Ref. 2]. 



Figure 2-1. Subsystems of the PHALANX Close-in Weapon System [Ref. 2]. 


Below this radar is the heart of the weapon system. It is the six-barrel 
M61A1 gatling machine gun as shown in Figure 2-2. It fires 20 mm (0.79 in) 
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tungsten (formerly depleted uranium) penetrators at a velocity of 1100 m/s (3600 
fVsec). An example of a projectile is shown in Figure 2-3. The penetrators are 
designed to mostly maintain their shapw as they hit and pass through the target. 
In this way, the penetrators can transfer their high momentum to the target in a 
concentrated manner, rather than fragment and perhaps partially reflect off the 
target. In experimental studies conducted at the Naval Surface Warfare Center in 
Dahlgren, Virginia, it takes on average six hits to destroy an incoming missile. 
Tests are currently being conducted against patrol boats, but no data is available 
yet. 



Figure 2-2. M61A1 PHALANX gatling machine gun [Ref 3.]. 
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Figure 2-3. Mk 149 PHALANX round [Ref. 2]. 

The gun has a selectable firing rate of either 3000 shots per minute or 4500 
shots per minute. A typical missile engagement can last from three to five seconds 
so that the higher firing rate is used more often. This means that somewhere 
between 225 and 375 bullets are fired at the missile. For engagements longer than 
five seconds, the barrels begin to deteriorate significantly due to the high barrel 
temperature. Therefore, this type of firing burst is strongly discouraged. The 
ammunition magazine is suspended below the gun and has been upgraded to hold 
1500 rounds allowing for roughly three engagements. Only one barrel fires at a 
time near the top of the barrel group while the spent casing extraction and breech 
loading occurs during the rest of the rotation cycle. 
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The effective range of the gun is about 2000 m (2187 yds). However, the gun, 
after calculating the target speed, begins firing when the target is about 5 seconds 
&om the firing cutoff point at 100 m (328 ft) from the ship. If the target is moving 
fast, the gun will fire so that the rounds will begin intercepting the target at 2000 
m. If the target is moving slow, the gim will wait until the target is 5 seconds from 
the ship even though this range could be much less than the maximum 2000 m 
range. The firing algorithm will be different for the anti-patrol boat surface mode. 

Beneath the gun is a cubical structure which houses the computer, radar 
electronics, and associated cooling systems. This is the brain of the system. It 
determines the threat, time to impact, where to point the gun, when to start firing, 
how to move the gun based on observing both the target and outgoing rounds, and 
when to cease the engagement. It is a closed-loop spotting system in that the 
computer receives feedback on how close the roiinds are coming to the target and 
calculates the necessary gun movements to ensure intercept. Figure 2-4 illustrates 
the concept of close-loop spotting. The fire control computer calculates aim error by 
measuring the projectile miss angle at the target and subsequently corrects the 
gun’s aimpoint. When the fire control computer determines that the target is 
destroyed, that is, when the target’s speed or radar cross-section drops below a 
certain threshold, then the firing stops and the gun slews back to its horizontal 
position and resumes the search mode. For ships with two or more mounts, they 
are able to talk to each other in order to prevent engaging the same target. [Ref. 
2 ] 

PHALANX defines an air target as something heading towards the ship, or 
something that could be heading towards the ship after a 9 G turn. This precludes 
the possibility that the gun will shoot down the ship’s own outbound missiles. The 
speed of the target has to be above a certain limit, so that firiendly helicopters, and 
such, are not accidentally engaged. The PHALANX gun is designed to protect the 
ship on which it is mounted, so it will not attempt to defend other ships. 
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Figure 2*4. Geometry of PHALANX missile engagement [Ref. 4]. 

B. PHALANX ROUND DISPERSION 

Based on dispersion data provided during a data gathering trip to Hughes 
Missile Systems Division, Pomona, California, it was learned that the PHALANX 
weapon system knows where the target is, but has had a hard time hitting it 
effectively as evidenced by the dispersion pattern shown in Figure 2-5. Here, one 
can see the six individual barrel dispersion patterns, which are generally good 
individually. However, the target is the object at the center of this particular 
dispersion pattern where virtually no rounds have fallen. The fire control system 
is doing its job by keeping the target at the center of the gun’s dispersion pattern. 
It appears that much of the dispers'‘on can be associated with the six barrel-gun 
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system. To elaborate on dispersion. Figure 2-6 shows the geometry used to define 
dispersion. The probability of a penetrator passing through a circle of radius r and 
standard deviation o is given by the distribution, 


p(r) = 


- 2o* 
e 

2na^ 


( 2 . 1 ) 


It is assumed that all the penetrators pass through a circle of infinite radius, so 
that, 

2ii 

rp(r)dnci8 = 1 (2.2) 

The probability that a penetrator will pass within a circle of radius one o, P^g, can 
be determined by integrating equation (2.1) in polar coordinates from r =0 to r = o. 


P<a = l^pirydrdB 



^0 2ao* 


rdr2n 


(2.3) 


Simplifying and solving the equation yields, 



= + 1 = 0.39 . 


Therefore, 39 percent of the penetrators fall withing one a while 61 percent fall 
outside. In the dispersion example of Figure 2-6, 39 percent of the penetrators of a 
PHALANX gun will pass within a circle of one o = 1.2 milliradians. Part of the 
dispersion ('1 milliradian) can be attributed to spin induced yaw after the biillet 
leaves the barrel. But, a significant part of the dispersion is directly attributed to 
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barrel tip motion. The observed 2.4 milliradians of dispersion corresponds to 3.6 
mm of barrel tip displacement of a 1.52 m (59.8 in) long barrel. 



.TOOPBeSENTEDAHEAAT AZIMUTH (MRAO> 

(000 FT WITH 200 FT OFF RANGE 



Figure 2*5. Experimental dispersion pattern of the PHALANX gun. Note 

that the target drone is at the center while the individual barrel 
dispersion patterns fall all around it [Ref 3]. 

Other dispersion patterns shown in Figure 2*7 show a less pronounced 
doughnut-shaped pattern than in Figure 2-5, but it is still generally visible, with 
the target remaining unscathed in the center. The dispersion patterns were 
measured using an (Dehler and battenboard system at a test facility in China Lake 
as shown in Figure 2-8. An Oehler system consists of a set of four microphones 
which will listen for the acoustic shock wave of the parsing supersonic projectile. 
The arrival time of each microphone is compared with the other microphones and a 
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39% of rounds in burst pass thru this circle 


Figure 2-6. Illustration of radial dispersion as defined for the PHALA>TX gun 
[Ref. 3]. 

bullet position is calculated. The battenboard follows the Oehler system and is 
used to record the roimd pattern and to distinguish between the pusher plugs, 
which drive the ro\ind through the barrel (see Figure 2-3), and the rounds 
themselves. This observation set-up was designed to test the gun dispersion and 
not that of the ammunition. 

The PHALANX dispersion pattern measurements also show that the six- 
barrel system gives a dispersion pattern greater than that of a single barrel fired 
separately. Thus, the structural dynamics of the six-barrel system appears to 
warrant an in-depth analysis. According to the Hu|^es data, the PHALANX gun 
system has a typical dispersion range of 1.2 to 1.5 milliradians and a maximum of 




Figure 2-7. A less pronounced circular PHALANX dispersion pattern [Ref 3]. 

1.6 milliradians. The U. S. Navy would like to achieve an average dispersion of 
less than 1.2 milliradians. In order to achieve this goal, the dynamics of the gun 
system will have to be better understood before any type of fix can be fully 
implemented. This is the premise behind this study. What effect the observed 
dispersion has on target destruction is the subject of the next section and was the 
motivation for imdertaking this project. 
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Figure 2-8. Experimental set-up for observing the PHALANX dispersion 
patterns [Ref. 3]. 

C. THE BALLISTIC DEBRIS PROBLEM 

Although the maximum effective range of the PHALANX g\m is 2000 m, field 
and computer simulations have shown that the destruction of an incoming missile 
does not occur, if at all, until after the missile is within 1(X)0 m (3281 ft) of the 
ship. On the surface this does not appear to be a problem. As part of a class 
computer simulation project, the resulting debris firom the destroyed missile 
continues on a ballistic tr^ectory toward the ship at roughly the speed of the 
missile. The firagment sizes, which can be on the order of a few to a hundred 
pounds, coupled with their high velocity, which is typically from 1000 m/s to 1500 
m/s, would cause a great deal of damage to an unarmored superstructure which 
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usually houses sensitive electronics, antennas and not to mention sailors. Metallic 
particles at 1000 m/s would easily punch holes through several bulkheads before 
being stopped. Although the PHALANX gun "worked" in destrojdng the missile, 
the ship would still suffer tremendous damage. A solution to this problem, would 
be to obviously destroy the missile at a greater distance from the ship. In order to 
achieve this goal the dispersion pattern would have to be tightened so more hits 
occurred at near the beginning of the PHALANX engagement rather than at the 
end. 

To clarify these points further, a computer simulation is used to describe 
fragment trajectories. In this model, a missile is assumed to be destroyed after one 
or several of the PHALANX rounds detonates the warhead of the missile. The 
range at which this detonation occurs is varied so that the number of fragments 
that hit the ship can be determined. The fragment size, shape and mass is 
approximated by distributing the fragment terminal velocities by the exponential 
distribution. 


-»r 



where is a fragment terminal velocity, and is the characteristic terminal 
velocity which is derived frrom typical fragment attributes. The terminal velocity is 
calculated by the expression 



(2.5) 


where m is the mass of the fragment, g is the acceleration of gravity, C is the 
coefEcient of air drag, p is the density of air and A is the debris fragment surface 
area. Additional velocity components caused by the force of the explosion are 
assigned to the fragments along each of the three coordinate axes. These velocities 
are determined from a Gaussian distribution with a standard deviation of 0 ^, in 
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each direction. The value of = 10 m/s was provided by the Naval Surface 
Warfare Center, Dahlgren. When the values of the variables described above are 
entered into the simulation program, the results appear like that of Figure 2-9 
This figure shows a typical best case scenario where the missile blew up at the 
maximum PHALANX range of 2000 m. It is divided up into a side view, top view, 
and an end view showing how the missile disintegrates as viewed from the ship. 
The missile is destroyed at the left most edge in the top and side views. The 
"ship" is represented as the right hand edge of these views. The black dots are the 
missile debris fragments and the black lines represent their trajectories. The 
"ocean" is at the bottom edge of each view. It is easy to record ship hits as those 
that reach the right hand edge. In this particular case, the missile was traveling 
at a fast speed of = 1000 m/s (3281 fl/s) and a height of 10 m above the water 
when it exploded. NF stands for the number of fragments. Here 6 missile 
firagments hit the ship with a velocity of around 100 m/s (328 fVs). 

In Figxire 2-10, the missile has the same characteristics, but now does not 
explode until it is intercepted by the last PHALANX bullet at a range R = 500 m 
(1641 ft) from the ship. The resulting debris hits the ship in a "shotgun blast" 
pattern where 18 of the 20 fragments reach the ship. This latter case is the more 
likely of the two scenarios for a nonmaneuvering, medium speed missile [Ref 5]. 

If the missile was more sophisticated in that it "jinked," or performed random 
maneuvers, or if it had a higher speed, then the intercept range is even more 
reduced. Therefore, it is readily apparent that significant debris from an effectively 
engaged missile is going to hit the defending ship in virtually every case. Such 
debris would certainly not destroy the ship, but might seriously impair her ability 
to fight. 

In other simulations, performed by Professor W. B. Colson, at the Naval 
Postgraduate School, the effectiveness of the gun in destroying an incoming missile 
is tested. In this simiilation, the round dispersion is assumed to have a Gaussian 
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Figure 2-9. Twoectories of fragments of a missile destroyed by a PHALANX 
gun at a range R = 2CX)0 m. Notice 6 of the 20 missile fragments 
hit the ship which is at the right hand edge of each view. 

distribution with a standard deviation using a typical PHALANX dispersion value 
of 50 = 2 milliradians. Figure 2-11 shows the results of a PHALANX engagement 
with an incoming missile which has a speed of = 600 m/s (1969 fVs) and a 
physical cross-section of 0.20 m (7.9 in). In thiR scenario, the missile speed is not 
unusual, but it does begin a terminal 5 G maneuver at a range of roughly 2000 m 
(6562 ft). The PHALANX gun destroys the missile at a range of roughly 200 m 
(656 ft) after it was only hit three times. The number of hits required to destroy 
the missile is chosen randomly and is better than the usual six. The effect of the 
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** Oabrla Slaulatlon: Part III ** 

V BlOOOn/s NFs20 C«1 psl R«500m 

n 



Figure 2-10. An incoming missile is destroyed by the PHALANX gun at a 
range of 500m. Most of its debris hits the ship. 

maneuver was to allow the missile to penetrate to within 200 m of the ship when it 
exploded and sends almost all of its debris at the ship. 

Another engagement is shown in Figure 2*12 where the missile does not have 
a terminal maneuver, but has a higher speed of 800 m/s (2625 ft/s). In this case, 
the missile is also destroyed, but like before, it is not destroyed far enough away 
from the ship so all of its debris strikes the ship. [Ref 5] 
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Figure 2-11. A terminal missile maneuver defeats the PHALANX gun [Ref. 5], 

These studies provided the motivation for trying to improve the dispersion of 
the PHALANX gun. In order to make such improvements, the gun dynamics have 
to be first understood. This study attempts to provide a solid foimdation upon 
which to study the gun. In an era of reduced defense spending, but not necessarily 
one of reduced potential for hostilities, the PHALANX gun is going to remain a 
ship’s last layer of defense into the next century. 
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Figure 2-12. A high speed, nonmaneuvering missile defeats the PHALANX 
gun [Ref. 5]. 
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III. FINITE ELEMENT THEORY 


A. THE DECISION TO USE FINITE ELEMENT MODELING 

In order to investigate the mechanical dynamics behind the dispersion of the 
PHALANX six-barrel system and to study the various design factors which 
influence projectile trajectory and accuracy, it was decided that a finite element 
computer model would be the best method for gaining insi^^t into this complex 
system. Many aspects of the gun dynamics would have been difficult if not 
impossible to understand by studying the gun itself during its operation due to its 
hiyh rate of fire and extremely high barrel temperatures as well as with problems 
of on-sight availability. Even high speed film does not provide as much benefit or 
nearly the flexibility as a computer model. A finite element analysis can provide a 
thorough model of the gun where various driving forces can be applied and the 
gun’s response, such as barrel tip motion, calculated. This also allows individual 
components to be studied in detail to understand their affect on the overall 
dispersion pattern. Furthermore, this finite element computer model is adaptable 
enough to support redesigning the PHALANX barrels and perhaps the thrust 
bearing, as well as further modifications, such as a change in projectile size or the 
addition of a gun barrel restraint. Finally, a computer code provides an economical 
solution by leapfi:x)gging past design changes made possible by computer simulation 
instead of expensive hardware tests. 

The finite element computer code called I-DEAS (Integrated Design 
Engineering Analysis) developed by SDRC (Structural D}mamics Research 
Corporation) was chosen to model the PHALANX gun due to its existing on-site 
license with the Naval Postgraduate School and its "completeness." That is, it 
precluded the need for other software packages allowing for a much smoother 
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transition between model construction and testing. I*DEAS is made up of the 
following set of modules: Solid Modeling, Drafting. Finite Element Modeling and 
Analysis, System Dynamics, Test Data Analysis, and Manufactiuing. This study 
uses a two-step approach. The first step uses the Finite Element Modeling and 
Analysis Module to develop a finite element model of the gun. This module allows 
normal modes, stresses, and static defiections to be calculated after the structure’s 
material properties and loads are defined. In the second step, the completed finite 
element model is transferred to the System Dynamics Analysis module where an 
in-depth forced response analysis is then performed. 

The primary purpose behind the finite element modeling of the PHALANX 
gun is to better understand some of the causes behind the apparent dispersion 
pattern. Specifically, barrel tip displacement pla3rs a critical role in the dispersion. 
One of the best ways to predict the motion of a system while being excited is to 
first solve for its normal modes. Normal modes of vibration have the iiseful 
property of being able to describe any complex motion (due to arbitrary initial 
conditions) as a superposition of normal mode amplitudes. So what are normal 
modes of vibration? Briefly, every structure, with its inherent mass and stiflhess, 
has its own distinctive set of normal modes. A structure vibrating at one of its 
natural fii^quendes does so with a characteristic pattern of amphtude distribution 
called a normal mode of vibration. A more mathematical explanation will be 
provided as the theory behind finite elements unfolds in this chapter. After 
modehng a structure as a series of finite elements and connecting points, or nodes, 
the computer code solves for these normal modes of vibration. Then, this model 
with its normal modes is ireitiferred to the Systems Dynamics Analysis module 
where excitation forces and damping can be applied. The response, such as barrel 
tip displacement, can then be calculated. A fourier transform of the time response, 
which performs a transient integration, estimates which normal modes are excited 
and their relative degrees of excitation. This same information can be revealed 
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using a frequency response function which is the method presented in Chapter VI. 
This latter method has the added benefit of not requiring a fourier transform. The 
resulting dispersion could then be modified, for example, by tuning these normal 
modes to different frequencies. 

The finite element method provides a discrete representation of a dynamical 
system [Ref. 6]. It involves the systematic solution of partial differential equations 
allowing the structural analysis of the complex systems to be easily performed by 
digital computers. Developed in the late 1950’s, it has been made more accessible 
with the advent of powerful desktop computers. It is especially iiseful for 
structures that contain a large number of degrees of fii^dom, where a variety of 
materials are used or where discontinuities such as thickness variations, holes, or 
curved beams are present [Ref 7]. These types of structures make exact solutions 
virtually impossible. 

The finite element method breaks down a structure into simple, discrete 
pieces such as beams, plates, or bricks, each of which is responsible for carrying 
some portion of the load applied to the structure. The equations describing the 
motion or elasticity of these discrete elements can be solved individually with 
boundary conditions matched at their nodal points to get the overall dynamic 
response of the structure. The calculated solution is not exact, but, with good 
modeling practice the error can be reduced to tolerable levels. Moreover, there is 
no other way to solve the problem. For displacement based finite element 
modeling, the error decreases monotonically. The finite elements and nodes are 
assembled to correspond as close as possible to the geometry of the actual 
structure being modeled. More elements are assembled in areas where the 
greatest stress gradient or curvature is expected or where the greatest accuracy of 
solution is needed. This often leads to a trade-off between computational speed 
and computer storage on the one hand and the necessary accuracy on the other. 
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B. THE DISCRETIZATION OF A SIMPLE ROD 


As an illustration of the finite element method, consider a single rod element 
of length L as shown in Figure 3.1 below. It has a node at each end point whose 
coordinates in the element’s local coordinate system (x ^) are x = 0 for node 1 and 
X sL for node 2. This element is situated in a global coordinate system (x,z) 
which would be used to describe a structure of elements. 



I I-1 I 


Figiire 3-1. A single rod element. 

To structurally represent the rod’s continuous distribution of mass, it will be 
broken down into an infini te number of discretized mass units dm. These mass 
units are then idealized by tracking their positions to the coordinates at the end of 
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the rod, specifically its nodes. Finally, an axial force is applied to the element 
along the rod’s local x • axis which will cause an axial deformation. 

In order to investigate the forces and displacements of this rod element it will 
be beneficial to use a displacement field variable ( u =u{.x^)) to describe the 
positions of the nodes at any instant in time. This is important since it is at these 
node points that the rod element is coupled to the rest of the structure. Any load 
information transmitted firom the rest of the structure is applied to the rod element 
at its nodes. Therefore, the response anywhere along the rod can be described in 
terms of the response at the element’s endpoints, that is, its nodal degrees of 
fi'eedom, which are finite. This is a finite discretization of an infinite element 
system. An exact solution of the displacement field is difficult considering that the 
force and moment distribution as well as the botindary conditions at the ends of 
the rod are not known. However, an approximation can be made by considering a 
finite number of degrees of fi:eedom for the displacement fields of each element 
[Ref. 7]. The coordinates for nodes 1 and 2 at any instant in time are respectively, 

«i«u(0,/) (3.1) 

and 

U2-u{Ljt) . (3.2) 

Since the load is an axial force, only displacements along the loc*! x • axis are 
considered, so the rod deformation can be described by these two displacement 
degrees of fineedom, u x and u 2 - 

The displacement field is a linear motion in x since there are two axial 
degrees of fii^edom [Ref. 7]. It is written as 

u(x,t) = Cl + C 2 X , (3.3) 

where Ci and C 2 are vinknown constants. Substituting equation (3.3) into equations 
(3.1) and (3.2) yields 
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«l = Ci 


(3.4) 


U2 = Ci + CaL . 

which forms a system of simultaneous equations 


1 0 


Cl 


ui 

1 L 




.“2 


where 


1 0 

-1 

1 

L 0 

1 L 

” L 

-1 1 


Solving for the Cj coefficients gives 


Ci = Ui 


and 


C2 = 


(—tt 1 + ^ 2 ) 

L 


Substituting equations (3.8) and (3.9) into equation (3.3) 


u(x^) = Ui 


uix 

L L 




X 

‘-T 

^ J 

Ui + 

L 

» 4 


= Oi(x )u i«) + a^x )u ) 


where 

and 


(3.5) 


(3.6) 


(3.7) 


(3.8) 


(3.9) 


(3.10) 


(3.11) 


02(x . 


(3.12) 


or in matrix notation 
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where 


r ^ 


u (x ,< ) = 


Ctl 


“1 

“2 


(3.13) 


a 


•^ttiCx) ajCx)!' 


(3.14) 


The response anywhere along the rod can now he described in terms of the nodal 
coordinates and the interpolation functions or shape functions ot^(x). These 
functions interpolate the element displacement between the nodes according to the 
boundary conditions and movement of the nodes. They also form a basis set where 
any motion of the rod can be described as a linear combination of these shape 
functions. 

The boimdary conditions defined earlier (see Figure 3-1) for the element 
impose the following requirements on the shape functions, Oj: 

ai(0) = 1 aa(0) * 0 (3.15) 

ai(L) = 0 aa(L)=l , (3.16) 

where each shape function represents the element shape when a unit displacement 
is applied. Looking at Figure 3*2, it is apparent that the shape functions can now 
be written explicitly as 

ai(x) = 1-L 


a^x ) * ^ 

The shape fimctions do not give an exact solution of the displacement field 
between the nodes, but as the niunber of elements are increased along the length 
of the rod a better solution is obtained. Referring to Figure 3-3, an exact solution 
of the displacement field can be written as 
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a^L ) = 1 


ai(x) = 1 - L 

ai(0) = 1 



0 X 





Figure 3-2. Shape functions . 



(3.17) 


where x’ and y’ are the coordinates of a point P along the length of the rod and the 


quantities 


and 


are the basis vectors. Therefore, any point along the length 


of the continuous rod can be described in terms of the basis vectors. However, it is 
the quantities x’ and y’ which are unknown and must be determined according to 
the motion of its endpoints. This is a linear assumption that is only valid if the 
transverse displacement is small, that is, x/Lcl. 
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Figure 3-3. Diagram of a point P along a deformed rod. 

C. MATRIX FORM OF EQUILIBRIUM EQUATIONS 

Since each normal mode of vibration is made up of the periodic exchange of 
kinetic and potential energies within the structure, it is necessary to determine the 
kinetic and potential energies for each element [Ref. 8]. Specifically, the mass and 
stiffiiess matrices of the element can be deduced after applying Lagrange’s 
equation to the internal kinetic and potential energies. In the case of the rod, its 
kinetic energy T and strain energy U are respectively 

T = ^j^[u(x^)]^dx 


(3.18) 




where y is the density per unit length and E is Young’s moduliis of elasticity and it 
is assumed that the rod is uniform with a constant cross-sectional area A. The 
strain energy occurs due to any axial elongation of the rod where du(x ,t) i dx is its 
corresponding strain rate of change. 

Differentiating equation (3.10) with respect to time yields 

u = ai(x )u ) + a2(x)ii2«) (3.20) 

and with respect to z gives 

4^ = a'i(x)ui(0 + a'2(i)n2(n . (3.21) 

dx 

Now substituting these expressions back into equations (3.18) and (3.19) defines 
the rod’s internal kinetic and potential energies respectively as 

V 

T = (oi(x )u 1 + 02(x liiate (3.22) 

t; = M J^V'i(x )u 1 + aVx )u ifdx . (3.23) 

The equations of motion can be derived firom T and U by applying Lagrange’s 
equation. It is defined as 

= 0 , (i = 1,2,3,...n) . (3.24) 

dt dq. dQi 

where L is the Lagrangian of the system. For conservative systems, where the 
sum of potential and kinetic energies remains constant, and mechanisms of energy 
loss such as internal firiction are ignored, the Lagrangian is the difference between 
kinetic and potential energies: 

L = T -U . (3.25) 

It is important to note that damping as a form of energy loss can be added into the 
equations of motion later (see Chapter VI). The generalized coordinates q. 
(i=l,2,3,...n) describe the configuration of the system, with n degrees of freedom 


31 




and hence n equations [Kef. 9]. In the case of the rod element, the generalized 
coordinates are 


= Ui Q2 - ^2 > 


so that 


<7 


.“2 


(3.26) 


Now applying Lagrange’s equation 

[ai(i )u 1 + 02(1 )ii 2 ]oii(i )dx (3.27) 

= EAj [a'i(x )u 1 + aVx )u 2 ) 0 'i(z )dx (3.28) 

au \ 

Therefore, the element mass matrix coefficients can be written as 

niij = yj^ Oi(x )aj(x )dx , (3.29) 

and the stiffiaess coefficients matrix as 

kij =EAj^a'i(x)a'j(x)dx . (3.30) 

So the equations of motion can be written in the familiar Newtonian mechanical 
forms 


>l 


*11 

*12 

Ml 

F2 

• 

*21 

*22 

.“2 






• ■ 

^1 

.-^2 

= 

mil mi2 

m2i ^22 


Ml 

“2 







(3.31) 


(3.32) 


It is important to note that equations (3.29) and (3.30) have the general 
quadratic form with n variables qi of the type 




(3.33) 





where the aij (i.e. y and EA) are coeffidenta independent of [Ref. 9]. This form 
is diagonal if it contains terms with squares of Qi namely anq;, 022^2 > > 

no cross terms like a i 2<719 2 9]- 

However, if there are cross terms in either the mass or stiilhess matrices, 
then there are a set of n coupled differential equations which describe the motion 
of the system. In the physical sense, this means that the motion of one node 
influences the motion of the other through these coupling terms. They are linearly 
independent. So the problem now is to simultaneously diagonalize the mass and 
stifBaess matrices if they are not already diagonalized. It is possible to exdte a 
normal mode of one node withou , exdting any other modes. The set of vectors that 
diagonalize the modal matrix are called "normal modes." This diagonalization is 
achieved by performing a linear transformation involving the orthogonal modal 
matrix 0, whose columns are the eigenvectors of the matrix K with respect to M of 
the eigenproblem 

iK - fl)2Af)0 = 0 . (3.34) 

Once the matrices are diagonalized, and hence uncoupled, then the n equations of 
motion are easily solved for less than n, rather than solving for n coupled 
(simultaneous) differential equations, and more importantly, each uncoupled 
differential equation is one of the normal modes. This process is called a transition 
to normal coordinates. It is important to note that the physical nature of the 
problem has not been changed, rather, the normal coordinates merely capture the 
physical configuration of the system so as to greatly simplify the mathematical 
description of the system. 
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D. FREE VIBRATION IN A STRUCTURE OF ELEMENTS 


Now that the equations of motion of a single rod element are understood, a 
whole structure can now be composed of elements like these. The equation of 
motion for an assembly of element matrices is 


[Mix +[iini =0 


(3.35) 


where [M] and [K] are square matrices and the quantities with a single underline 
imply a column matrix. As an example, if the original continuous rod is 
represented by eight connected elements instead of one, as shown in Figure 3-4, 
each one with its own 2x2 stiffness matrix, Ki , then the eqxiation of motion would 
appear as 




\k,] . . . ■ 

“1 

F 


. [^^23 . . 

“2 

L J 


• * • • 

• 



. . [Xsl 



I 

(3.36) 


where each element stif&ess matrix, Ki, has an overlapping stifiness with adjacent 
elements at their common degrees of freedom. 

The general method for solving the free vibration problem, that is, finding the 
eigensolution of (3.35) is to use a solution of the form 


x«) = Cie*“' 


where x(t) is the column matrix 


xit) = 


Xiit) 

x^t) 


(3.37) 


(3.38) 


describing the position of the element end points or nodes. The solution consists of 
a vector spatial part. 


$ = 


♦2 


(3.39) 
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Figure 3*4. Continuous rod divided into 8 elements. 


scalar temporal part, e‘“* and a complex constant C. This form of solution was 
used since the mass and stiffness associated with each element is similar to a 
mass/spring system which executes simple harmonic motion. 

Taking the necessary derivatives of the solution and substituting back into 
the equation of motion yields the homogeneous equation 

(-a)2[Af 1$ + [K}^)C€*^ = 0 . (3.40) 

If C = 0 then there would be no motion, so that possible solutions to this equation 
have to come firom 


([if]-©^[Af 1)^ = 0 , 


(3.41) 





which is an eigenvalue problem. Part of the process in finding solutions to 
eigenvalue problems is to take its determinant 

det([K]-<aHM])^0 , (3.42) 

which is called the "characteristic equation." This is the necessary and sufficient 
condition for these n homogeneous equations in n unknowns to possess a non¬ 
trivial solution. The values of which render ([/f] - (o*[Af])"^ singular are called 
eigenvalues. For two degrees of freedom, there will be two roots of this equation, 
coi and (Os. The nodes can vibrate in simple harmonic motion with only these two 
frequencies. To determine the actual motion of the system, Oi and (O 2 are 
substituted back into equation (3.37) to find Xi(t) and X 2 (t) for each frequency. 
Since the equation of motion is homogeneous, the coefficients Ci and C 2 are 
undetermined which means the system fan vibrate with an arbitrary amplitude for 
either characteristic fi?equency. The relation between and C 2 is essentially a 
phase relation. [Ref. 9] 

The next step is to apply this method to the whole system of elements with n 
degrees of fireedom. Starting with the eigenproblem defined earlier 

([fir]-(i>2[Af])$ = 0 , (3.41) 

except now [K] and [M] are made up of a whole series of two by two matrices, one 
for each element. For systems with many degrees of fireedom and especially with 
cross terms appearing in the mass and stifEhess matrices, the n coupled 
differential equations describing the motion now become difficult to solve. 
However, a relatively easy method to solve such a complex system is to 
simultaneously diagonalize the matrices [M] and [ff]. Then, it is trivial to find the 
eigenvalues and eigenvectors. In this case, the diagonal elements are now the 
eigenvalues and the eigenvectors are then the columns of the matrix used to 
diagonalize [M] and [iT]. 
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The method used to render [Af] and {K] diagonal is a transformation to 
normal coordinates. The problem is to find an orthogonal matrix* d> which 
simultaneously diagonalizes [Af ] and [AT], If [<I>] is orthogonal then for unity modal 
mass normalization 

[<S>f[M]m^[I] , (3.43) 


or 


[<i>]^[Af] = [d>r^ 

Recall that the identity matrix can be defined as 

0 *^ 4 ^ * 0^0 “ 2 


where 

10 .. 
0 10 . 
^ . 0 1 0 
..01 


(3.44) 


(3.45) 


The process of simultaneously diagonalizing the [Af] and [K] matrices can 
begin by rearranging the eigenproblem equation (3.35) as 

[An$ = a)*[Af]$ . (3.46) 

Corresponding to the ith and jth normal modes are the vector equations 

[K]^^ * Wj^lAf , (3.47) 

and 

[K]^ = . (3.48) 

Using the identity matrix to make them scalar equations, equation (3.47) is 

•T 

multiplied by ^ , where the transpose changes the n x 1 column matrix into a 
1 X n row matrix, and equation (3.48) is multiplied by so that 

=a)iV’'[Af]<^‘ , (3.49) 
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and 


. (3.50) 

Since [K] and [M] are real, symmetric matrices [X"] = [A’]^ and [Af] = [M]^. 
Making use of this property, the transpose of equation (3.50) yields 

. (3.51) 

which when subtracted from equation (3.49) gives the scalar equation 

0 = ((i)i*-a)/)^’‘lM)$‘ . (3.52) 


If (Oi , there are two different normal modes, then 


♦i W 


^11 ^21 
A/21 ^^22 


♦ l 

♦2 


= 0 


(3.53) 


One of the key characteristics of normal modes is that every normal mode is 
orthogonal to every other mode. Specifically, 


r 



i *j 
i • 


(3.54) 


The mass matrix [M] is now diagonalized into a new modal mass matrix M^. This 
matrix M} is a measure of how much mass is participating in the ith normal mode. 
Similarly for the non-diagonal stiffiiess matrix [K] 


= 


0 i 
K, i =J 


(3.55) 


Therefore, diagonalization is achieved by a linear transformation using the 
orthogonal modal matrix <1>, now defined as 


[<!>] = 


W] (♦' Jf 


(3.56) 


which represents the spatial and temporal parts of the mode shape. It is composed 
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of the eigenvectors and V. The eigenvalues form a spectral matrix 


0 

0 ( 0 * 


Returning to the equation of motion 

[Af]x =0 . (3.35) 

and introducing the orthogonal modal matrix as part of the solution 

X = [<!>]£ . (3.57) 

which when substituted back into equation (3.35) yields 

[MKOlj=0 . (3.58) 

To perform a linear transformation, this equation is multiplied by [<I>]^ so that 

[OfEJIfKOlsi=0 . (3.59) 

which diagonalizes [M] and [F\ to create a new normalized relation 


Mu 0 . . ■ 



I 

0 

0 M 22 0 

r • 


0 K 22 0 

0 . 0 

0 M(j 

[9 

■f 

0 . 0 

0 Ky 


= 0 


(3.60) 


Now the n coupled differential equations are uncoupled. Every ith normal mode 
now has its own ith row which can be lifted out, such as 


M|7i + Kj(7i = 0 

which is like the familiar harmonic oscillator mass/spring system where 


(3.61) 


or 


Qi + = 0 
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Any arbitrary motion of the system can thus be written as a linear combination of 
the eigenvectors 

= 9i[«i]+<?2[«2] . (3.62) 

where qi-Acosait and q 2 = Acos. Here A is an amplitude coefEcient based 
on the perturbing force. Hence, the finite element computer code solves for the 
normal modes of a structure by finding a modal matrix <I> which will diagonalize 
[ilf ] and [Kl 

Therefore a normal mode is made up of a particular eigenvalue (frequency) 
and a particular eigenvector (mode shape), which is orthogonal to every other 
normal mode. A normal mode describes how much mass and stif&iess is involved 
in the structure as it vibrates. That is, the coordinate is a measure of how much 
the mode shape d>, is present The value of finding the normal modes of a 
structure is that any oomplex motion can be represented as a linear combination of 
normal modes. For such an intricate system as the PHALANX gun, where there 
are many components contributing to its dynamic response, to find its normal 
modes suddenly makes it a far easier system to understand. That is, if the normal 
modes of the gun are obtained, which then forms a modal component model, when 
excitation forces are applied, the res\ilting displacements are the result of exciting 
certain normal modes. An analogy to this is that any color can be made up of 
various amounts of the primary colors red, blue and yellow. So in order to obtain a 
set of accurate normal modes the finite element model must be as close a 
representation of the real gun as possible which is the topic of the next chapter. 


•*1 

Jt2 
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IV. MODEL DESCRIPnON 


A. TYPES OF FINITE ELEMENTS USED IN THE MODEL 

At the outset of this project, the gun barrels alone were thought to be a 
primary source of dispersion. In fact, previous to this project, the Navy had 
already let a contract for a new barrel design. So, the initial proposal for research 
by the Naval Postgraduate School called for a finite element model of the six-barrel 
system. However, in consultation with Mike Hatch, a private-industry finite 
element expert, it was agreed that in order to truly represent the complex 
dynamics of the gim, it would be important to model not just the barrels, but the 
whole gun assembly through to its grounding points to the g\in mount. This 
holistic approach, while greatly expanding the scope of the work, would not only 
offer a much better assessment of the gun system, but would greatly expand its 
applications as well. The school would also be the first organization to model the 
complete gun assembly, thxis offering a complete model in one package. The 
benefits of such a model are many. Once the model was assembled and the 
dynamic response validated against the current dispersion pattern, new 
components could be "swapped out" in the model and their effect on dispersion 
would provide a valuable assessment. Not only could individual components such 
as the barrels be tested, but combinations of components could be tested together. 
This model therefore, would go a long way toward understanding the dynamics of 
the gun as well as providing a valuable tool for finding ways to tighten the 
dispersion pattern. 

Once blueprints of the model were received, the next big decision involved the 
appropriate selection of finite elements to represent each component. In all there 
are ten different parts of the gun, coimting duplications such as the barrels, as just 
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one component. The finite element modeling experience Professor Gordis and 
Mike Hatch was helpful in deciding which areas of the gun system should be 
modeled with great detail and areas where time-consuming attention could be 
safely passed over. Table 4-1 lists the components and the corresponding element 
choices. Before discussing the properties of the elements that were used in the 
model, it is important to give an overaU view on element selection. Each element 
comes with a price. Some, such as bricks, are more expensive to use than others 
since they require more nodal points to define them, in this case eight. Others, 
such as beam elements, as discussed in the previous chapter, require only two 
nodes. In order to shorten the computational time and make the resulting modal 
shape animation as smooth as possible, it is important to choose as many less 
expensive elements as possible. At the same time, to accurately represent the real 
gun, the expensive elements are sometimes the best choice. Hence, it is a balance 
that often goes through several iterations before a final choice is made. As can be 
seen in Table 4-1, the gun is comprised of five categories of elements, with many 
different variations within each category, such as changes in dimensions, material 
properties, etc. The following subsections describe the general properties of each 
category of element used in the model. 

1. Brick Elements 

As mentioned before, the solid brick elements are the most complex 
elements to use in that they require eight nodes to define. A single brick element 
is shown in Figure 4-1. A brick element is one of the most versatile elements to 
use since it allows modeling of any structure with any geometry [Ref 10]. The 
general geometry of the brick is that its width is roughly comparable to its length. 
It is a three-dimensional solid whose nodes are allowed to translate along each 
axis, but not rotate. Therefore, it has only three degrees of fireedom. It is best 
used to represent solid, invariant structures without many intricate details or 
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curvature, such as engine blocks or in our case the breech. Fvurthermore, more 
bricks are required to reflect bending and shear than do thin shells or beam 
elements. A structure represented with bricks often looks the most realistic. 
However, it comes with the price of reduced accuracy, especially if not enough 
bricks are used, and in any case, it is usually the most inefficient element to use 
[Ref. 11]. 


TABLE 4-1. ELEMENT SELECTION 


Elementary Model 


r—- 

Component 

Element Types 

Muzzle Clamps 

Bricks, Thin Shells, Rigid Elements 

Mid-Barrel Clamps 

Bricks, Rigid Elements 

Barrels 

Beams with tubiilar cross-section 

Stub Rotor 

Bricks, Thin Shells 

Rotor 

Bricks 

Needle Bearing 

Bricks, Springs 

Casing 

Thin Shells 

Thrust Bearing 

Bricks, Springs 

Recoil Adapters 

Bricks, Rigid Elements, Springs 

Ball Joint 

Rigid Elements, Springs 
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Figure 4*1. A solid linear brick element. 


2. Thin Shell Elements 

A thin shell or plate element is a two^limensional solid whose thickness 
is small compared to its length or in-plane dimensions. A thin shell element is 
shown in Figure 4-2. It is a much more efficient element to use than a solid brick 
since it only requires four nodes to define its shape. Its thickness is defined when 
the element is created, but it does not show up in the display. The best 
applications for it are in modeling situations where a wall or membrane has a 
nearly constant thickness or its thickness varies linearly from one edge to another. 
The element’s nodes are allowed to not only translate, but rotate as well, thus 
giving it six degrees of fireedom. These additional degrees of fireedom allow fewer 
thin shell elements to be used than brick elements in obtaining the same accuracy 
for bending and shear motion. One disadvantage in using thin shell elements is 
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that they are not able to model stresaes that vary through ita thickness as can 
brick elements [Ref. 11]. 



Figure 4*2. A thin shell element. 


3. Beam Elements 

Beam elements are the simplest elements in that they only require two 
nodes to define them. The cross section of the beam is defined when the beam is 
created and can have as detailed and complex a shape as the user desires. Figure 
4*3 shows a beam element as it appears in the display of the finite element model 
and its cross section which only appears when referencing the beam’s properties. 
This cross section is invariant along the element’s length. To represent a beam of 
varying cross-section, several elements are reqxiired. For example, the gun barrels 
are represented using beam elements, where the barrels’ varying thickness is 
reflected in the varying cross-section of 12 elements. The nodes are allowed the 
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full six degrees of fireedom which is usehil in displaying bending and deflection. 
Again referring to the modeled gun barrels, the beam elements provide a smooth 
anim ation picture which gives better insight into the dynamics of the gun. A 
possible problem with using beam elements is that they do not show stress through 
their cross section and the cross sections can be easily oriented the wrong way [Ref. 
11]. However, in most cases, there is more interest in the maximum, outer fiber 
stress than that of the inner stress. 



Figure 4-3. (a) A beam element, (b) A representative beam element’s 
cross section. 


4. Spring Elements 

Spring elements pan be defined with two nodes (node-to*node) or in terms 
of one node (node-to-ground). They can also be either a translational spring 
operating along a defined axis, or a rotational spring which represents the bending 
moment or torque that is required to rotate one node relative to another or to 
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ground [Ref. 12]. Figure 4-4 shows a node-to-node and a node-to-ground spring. 
The stifhiess is defined along each axis or un.^udally for a translational spring or 
about an axis for a rotating spring. Three rotational degrees of fireedom are 
assigned to each node. In the PHALANX g\in model they were used to represent 
the stiflhess of the ball bearings and the recoil shock absorbers. 


m-* 



Figure 4-4. (a) A node-to-node spring element, (b) A node-to-ground 
spring element. 


5. Rigid Elements 

Rigid elements are massless, infinitely stiif bars used to connect 
components together. That is, the rigid element implies that the motion of the 
nodes on the element are to be related to each other as if they were connected 
together by an infinitely rigid bar [Ref. 12]. Two nodes that are connected 
together with a rigid element move together according to a linear algebraic 
expression called a constraint relation. An example of such a relation is given 
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below. 


11-0.5x2 = 0 (4.1) 

In this example, the first term describes the distance node 1 moves in the x 
direction and the second term that of node 2 in the z direction. The right hand 
side of the equation, which is a constant term, must be zero when solving for the 
normal modes of vibration. This equation can be written in matrix notation as 


■Kl 


0.5 

*2 

.^2. 


1 

.*2 


Representing the coefficient matrix with T, this relation can then be written as 

^original = » (4.3) 

where the nodal position is defined in terms of a coefficient matrix and the 
constrained node matrix x«.. Substituting thin equation into the equation of motion 

F^KTz^ , (4.4) 

and multiplying by the transpose of the coefihcient matrix T yields 

TJF = r/XTxc , (4.5) 

which reduces to the constrained equation of motion, 

F^^K^x, . (4.6) 

This is an example of how constraint equations transform an unconstrained 
stiffiiess matrix A* to a constrained stiffness matrix Kg. A rigid element can be 
\ised to connect firom two to thirty-two nodes together. Figure 4-5 shows a rigid 
element connecting five nodes together. Each node at the element’s end points are 
allowed all translations and rotations that are allowed for that node. For example, 
if the rigid element is connected to a node of a brick element, then it is only 
allowed to translate. In the PHALANX gun model, rigid elements were used to 
connect the barrels to the clamps and stub rotor and the recoil adapter springs to 






the g\in body. "End release codes" can also be defined allowing the simulation of 
hinges or ball-and-socket joints [Ref. 12]. 



Figure 4-5. A rigid element connecting five nodes together. 

A further classification within each element category is the element order. 
That is, the number of nodes per element can be increased to increase the order of 
the equations used to interpolate the strain between them. There are two nodes 
required for each edge in a linear element, three for a parabolic and four nodes for 
a cubic [Ref. 11]. All the elements used in this model have linear order. 

B. CONSTRUCTION OF THE PHALANX GUN MODEL 

A breakdown of the components used in the construction of the PHALANX 
gun model is shown in Figure 4-6. This is a hidden line image which shows the 
geometry of the elements on the side of the gun parts facing the viewer. Figure 4-7 
is a shaded image of the completely assembled gun. Here, the individual elements 
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cannot be seen, but the effect of lighting and solid colors gives a realistic picture of 
the gun it represents. Both figures provide a brief overview of how the gun is 
assembled. The six barrels each slide through the stub rotor and are screwed into 
the rotor a distance of approximately 7.6 cm (3 inches). The stub rotor itself is 
bolted to the rotor and serves to support the barrels near their base. The six- 
barrels are connected together with a mid-barrel clamp and a forward muzzle 
clamp. This whole assembly is now attached to the gun body through a rear needle 
bearing and a double-roller contact bearing. The assembled gxm body is coupled to 
the rest of the CIWS mount at three locations. Two recoil adapters are positioned 
on either side of the gun body where they connect the gun body to the mount via a 
pivot joint. The third grounding point is a ball joint at the rear of the gun body. 
All of these components will be explained in further detail in the following 
subsections. 

1. Gun Barrels 

The six rifled barrels are approximately 150 cm (59 inches) long with an 
inner diameter of 2 cm (0.787 inches) and weigh about 8.2 kg (18 lb). The barrels 
are secured to the breech with interrupted locking lugs that mate after the barrels 
are screwed in past three threads [Ref. 13}. As mentioned before, the six barrels 
are represented using beam elements. Each barrel is made up of 39 elements with 
12 different cross-sections. Figure 4-8 shows a single barrel with its elements 
separated. 

2. Muzzle Clamp 

The muzzle clamp restricts barrel movement near the barrel tips. It 
consists of two plates held together with a steel web meshing. The clamp also 
restrains individual barrel movement during firing and a self-locking nut is used to 
secure the clamp to all the barrels [Ref. 13]. The plates are represented using 
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Figure 4-6. A hidden line image of the diaassembled components of the 
PHALANX gun. 

brick elements while the web meshing is modeled with thin shells as can be seen in 
Figure 4-9. 

3. Mid-Barrel Clamp 

The mid-barrel clamp is positioned roughly near the center of the barrel 
group. It physically locks the barrels to the clamp with tabs that fit into cutouts in 
each barrel. It prevents the barrels from becoming unscrewed in the rotor during 
firing. The mid-barrel clamp, as shown in Figure 4-10, is represented entirely with 
bricks. 
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Figure 4*7. A shaded image of the PHALANX gun. 


4. Stub Rotor 

The stub rotor which is bolted to the rotor supports the base end of the 
six barrels. It permits alignment of the barrels so that the locking lugs ran be 
engaged [Ref. 13]. like the muzzle clamp it consists of two plates separated by a 
metal meshing. Figure 4-11 shows the modeling of it with bricks and thin shells. 

5. Rotor 

The rotor is the most massive component in the gun body. It is 
hydraulically rotated at a speed of 750 rpm (12.5 Hz) in the coimterclockwise 
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Figure 4-8. A single barrel with with 39 beam elements. 

direction when viewed firom the rear [Ref. 13]. The gun assembly, which is shown 
in Figure 4-12, is bolted onto the front of it. The rotor is shown in FHgiire 4-13 and 
is made entirely of bricks. The rotor is coupled to the gun body through the rear 
needle bearing and the forward double roller contact bearing. 






Figure 4'9. The muzzle clamp represented with brick and thin shell 
elements. 

6. Double Row Contact Bearing 

The double row angular contact bearing, often referred to in this text as a 
thrust bearing, consists of two rows of ball bearings for a total of 76 balls. The 
inner race attaches to the rotor while its outer race attaches to the gun body. Its 
purpose is to align the rotor/gun assembly and transfer the firing load to the rest of 
the g\in body. It provides both radial and axial support. Figure 4-14 shows the 
races represented by bricks. The ball bearings are modeled using springs with the 
loaded contact angle being represented by an angle determined through an 
accepted computer program from Mike Hatch. 


54 









Figure 4-10. The mid-barrel clamp represented with bricks. 

After the decision was made to model the whole gun subsystem, careful 
attention was given to the blueprints of the gun body. It then became readily 
apparent that the two bearings were playing a critical role in the dj'namics of the 
gun. They provide the only means of support to the rotating and firing six-barrel 
gun assembly. Therefore, the intricate details of the bearings, especially the double 
row bearing, were modeled as can be seen firom a cutaway view showing the cross 
section of the double row bearing in Figure 4-15. It is also noteworthy that this 
bearing is mounted to the rotor through a narrow lip which was also accurately 
modeled and is pointed out in that figure. The entire firing and rotating forces are 
transmitted through this lip which seems to be subjecting it to a severe stress. 
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Figure 4-11. The stub rotor represented with bricks and thin shells. 

Obtaining the axial, radial and torsional stiilhess values for the double 
row contact bearing was difficult since the manufacturer, Kaydon Bearing in 
Muskegon, Michigan did not have records for this bearing in this application. 
However, after supplying Kurt Sheridan of Kaydon with the reaction force 
information, determined from the finite element model, he was able to provide 
approximate figures of the bearing stiffness and contact angle. With this 
information Mike Hatch was then able to provide spring stiffness values and an 
angle with which to mount the springs in relation to the inner and outer races. 
There are a total of 36 springs, one for each of the 18 elements used to define its 
circumference then multiplied by two to represent the double row of bearings. 
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Figure 4-12. The gun barrels with clamps and stub rotor. 

7. Needle Bearing 

The needle bearing provides the rear support for the rotor. It serves only 
as a radial support and not as a thrust absorber. The inner race of the bearing is 
attached to the rear wall of the gun body while the outer race is attached to the 
rotor. Figure 4-16 is a translucent picture showing the inner assembly of the gun 
body. In this figure it is possible to see both the needle bearing and the double row 
bearing and their relation to the rotor and gun body. Similar to the double row 
bearing, the bearing stiffiiess is represented with 18 springs. In order to correctly 
represent the actual bearing stifihess, ^6(aniv> ba effective spring stiffness, K^ff, 
must be calculated. No matter which way the rotor moves on the real gun, it 
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Figure 4-13. The rotor represented with bricks. 


experiences a constant bearing stiffiiess, whereas in the modeled gun which has 
only 18 springs, the angle between them must be taken into account. Figure 4-17 
demonstrates the geometry needed to make this calculation. Consider a force 
applied in any one direction to the model bearing. The equation of motion would 
be 


F^Ku^X 


(4.7) 


Here, Ktot is the total spring stifEaess value that would take into accoiint all 18 
springs. However, Ktu ^ Ff^taring since as it stands now, the modeled bearing will 
present a larger spring stiffness than the actual bearing. To correct for this 
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Figure 4-14. The double row contact bearing. 


problem, the 20P angle between the springy must be described by the relation 

Ktot = i,K,ffCoaai=Kb,ari^ . ( 4 . 8 ) 

>>o 


Solving for Kgff gives 

t-18 

JCgf f ^ COSOtj — '^bcarinf ' 
i -0 

where 


^bearing =Kgffi2c08(f + 4cO8 20® + 4c08A(f + . . .) 


Rearranging terms gives, 


( 4 . 9 ) 


( 4 . 10 ) 
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Figure 4>15. Cross sectional view of the double row contact bearing. 
Note the narrow lip which attaches to the rotor. 


= 


^bearing 

11.5175 


(4.11) 


Therefore, after obtaining the stifhaess value of the needle bearing from the 
manufacturer, it was divided by a factor of 11.5175 to fit the model [Ref. 14]. 


8. Gun Body 

The gun body is a steel housing which contains the bearings, rotor and 
stub rotor as shown in Figure 4*16. Since the wall thickness is only 0.5 cm (0.2 
inches) and is small compared to its other dimensions, it is represented by thin 
shells. It is part of the load path between the firing barrels and ground. On either 
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Figure 4-16. Translucent image of the gun body. It shows the relation 

between the two bearings and the rotor. Also visible are the 
recoil adapters and the rear support. 

side of the gun body, roughly opposite to the double row bearing, are the two recoil 
adapters. At the rear of the gun body is the third attachment to ground, the ball 
joint. 


9. Recoil Adapters 

The two recoil adapters serve to absorb a portion of the forces caused 
when the gun is fired and then the associated counter-recoil force [Ref 13]. Figure 
4-18 is a photograph of the PHALANX gun on the USS ARKANSAS which shows 
the construction and mounting of an actual recoil adapter. Since the overall shape 
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Figure 4*17. Spring placement for the needle bearing. 

of the recoil adapter plays a minor role in the dynamics of the gun, it was possible 
to represent its recoil spring stifhiess, placement and attachment to ground simply, 
but accurately. Looking at Figure 4*16, the recoil adapter’s mass and stiffness was 
modeled using brick elements in the form of "ears” sticking out on either side of 
the gxin body. In order to get the correct placement of the roller, rigid elements 
were xised to coimect the brick element ears to the recoil spring. In turn, the recoil 
spring was connected to a roller. The specifics of the roller will be discussed in the 
bovmdary conditions section. The stiffness value for the spring was obtained from 
a hysteresis curve as shown in Figure 4-19 [Ref. 15]. Here the slope for the upper 
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and lower branches of the curves was calculated and averaged to get the stiflhess 
value \ised in the model. 



Figure 4-18. Photograph of the gun body. Note the arrangement of the 
recoil adapters and their construction. 


10. Ball Joint 

The rear ball joint provides the third support of the tripod supported gun. 
The geometry of it is displayed in Figure 4*20. It consists of a 3.84 cm (1.5 inch) 
diameter ball supported by a 9.63 cm (3.79 inches) long rod 2.54 cm (1 inch) in 
diameter which in turn is attached to the gun mount. The ball fits into a socket 
mounted on the back of the rear wall of the gun body. When the gun fires, it sets 
back by about half an inch and oscillates about the new equilibrium point. The ball 
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Figure 4-19. Hysteresis cxirve of the recoil adapter stiflhess. 

slides into the socket which only provides radial support. Since this is a relatively 
simple arrangement, a separate side model was built to establish the stiilhess of 
the socket. This stiffiiess was represented by two springs along the y and z axes 
perpendicular to the length of the gun thus providing the radial support. The 
springs were attached to groimd on one end while their other end connected to 
rigid elements that attached to the gun body. These rigid elements represented 
the support given by the cup to the gun body. Figure 4-21 shows this arrangement 
in the computer model. 
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Figure 4>20. Blueprint of the rear ball joint. 

C. BOUNDARY CONDITIONS 

Having modeled the components of the gun, the next msyor step is to ensure 
that the proper boimdary conditions are reflected in the model. The boundary 
conditions define how the components are connected to each other as well as how 
the whole g\m is attached to the gun mount. The next category of boundary 
conditions that must be determined are the types of degrees of freedom. 

As the model stands now, it represents an idealized gun with no worn 
bearings, or recoil springs and no tolerance between connected components. It may 
be better than a brand new gun. As mentioned earlier, the barrels are attached to 
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Fissure 4*21. Representation of the ball joint in the model. 

the clamps, stub rotor and rotor with rigid elements. Figure 4*22 shows an 
example of a barrel connected to a clamp. The stub rotor is connected to the rotor 
with rigid elements as well The rotor is connected to the gun body through the 
two bearings. The gun body in turn is connected to ground via the recoil adapters 
and ball joint. 

Figure 4*16 expands further on the specific ranges of motion allowed at these 
three supporting points showing how they were defined in the model. At the two 
pins which hold the recoil adapters to ground in the real gun, they are allowed to 
rotate about the z or horizontal axis, but do not rotate or translate in any other 
direction. To properly represent this condition in the real gun, the node between 
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Figure 4-22. Example of rigid elements connecting a barrel to a clamp. 

the rigid elements and the springs is allowed to move in the z or axial direction 
and to rotate about the z axis. The translation in the x direction is necessary in 
order to allow the springs to be exercised. The springs in turn are attached to 
groimd in order to satisfy the real boxindary conditions, albeit broken up into two 
locations, given the unique representation of the recoil adapters in the model. 
Motion at the ball joint of the real gun is in the fore-all or x direction, with the cup 
providing restraint for translational motion in the y or z directions. All rotations 
are allowed, however due to the two forward supports this motion would be 
unlikely. In the computer model, the springs restrict motion in the y and z 
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directions, thus like the cup, and there are no restrictions for motion in the x 
direction. 

Unlike the real gun, this model does not spin. This is where kinematic and 
unknown degrees of freedom come into play. The normal mode of rotation about 
the X axis is prevented. Since the normal modes are being solved without 
gyroscopics, that is, without an unbalance (centrifugal) excitation, it is restricted in 
the model by using kinematic degrees of freedom. Four kinematic degrees of 
freedom were used to restrict the rotation of the rotor inside the gun body. This 
does not substantially degrade the understanding of the gun’s dynamic motion, 
since the excitation force occurs at the same place and is virtually independent of 
the rotation. That is, as the rotor is turning and the barrels are sequentially firing, 
they all fire at the same location. The barrels are not actually rotating in the 
computer model, but the excitation force is applied at the base of a barrel which is 
in the proper location to represent each barrel as it fires in the rotating gun. The 
unknown degrees of freedom correspond to the differential equations solved by the 
computer. In this model, there are 22,470 degrees of freedom. 

D. COMPARISON BETWEEN THE REAL AND MODELED GUN 

After completion of the PHALANX gun model, one of the first steps toward 
validation is to compare its properties with those of the real gun. Based on 
information provided by Dick Sirola, the properties used for comparison were mass, 
center of gravity, and moments of inertia for each of the three coordinate axes. 
This information and their percent difference between the computer model and the 
real gvm is presented in Table 4'2. Ftirther comparisons can be made and the 
model refined to provide a better match. These small differences are not important 
to the conclusions reached by this study. 
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TABLE 4-2. COMPARISON BETWEEN REAL AND MODELED GUN 


Comparison 


Mua (kg) Center of Gravity cm) Moment of Inertia (kg-m 




EM 


OH 

/= 

tyy 

— 

M61A1 

113.5 

53.83 

0.04 

0.0 

6796 

2.38e® 

2.38e® 

Model 

120.5 

55.2 

-0.304 

-0.053 

7240 

2.59e® 



HHIHH 







%Difif. 

5.81 

2.48 

5.75 

5.40 

6.13 

8.01 

8.40 


Another important area of oomparison is to determine whether the actual 
forces in the gun are linear or non-linear. The computer will only solve the linear 
problem. Experimental data verifies that the stiflhess calculations are linear for 
the range of barrel tip motion less than - 0.25 cm (0.1 inch) [Ref. 16]. Figure 4*23 
shows the test results of muzzle deflection versus a load up to 300 lb applied in the 
vertical direction at the barrel tip. Since the dispersion data shows that barrel tip 
displacement is less than = 0.25 cm, the system is considered linear. The linear 
computer solution is compatible with the experimental data. 
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Figure 4-23. Experimental data of muzzle deflection versus load. 
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V. MODES OF VIBRATION 


A. MODES OF VIBRATION OF THE GUN BARRELS 

One of the first components of the model to be developed were the 
barrels. Two different levels of modeling detail were used to create two different 
versions of the barrel. The first, called the advanced barrel, involved the use of 
496 brick elements and is shown in Figure 5-1. The brick elements were used to 
provide better detail in modeling the barrels. The second version, called the simple 
barrel, used 39 elements with 8 cross sections and is shown in Figure 5-2. The 
latter version, with its fewer elements, is the best choice for providing a smoother 
animation picture and keeping the number of degrees of freedom within a 
manageable budget. The advanced barrel version was kept for comparison 
purposes so that when future barrel designs are tested using this model of the gun, 
the most complete finite element modeling technique can be used. 

The first component of the PHALANX gun to be received at the Naval 
Postgraduate School was a single barrel. It provided an opportunity for a 
validation test of the gun model. The barrel was set up in a mechanical engineering 
lab, where it was suspended with two low-stiffiiess cords to determine its normal 
modes of vibration [Ref. 17]. This method is called the "firee-firee configuration." If 
the normal modes were found while the barrel was clamped at its base end, for 
example, the modes would change to match the clamped boxindary conditions and 
would not represent properties of the barrel above. In the firee-firee configuration 
the modes of vibration of the real gun barrel can be determined for use in the 
computer models. 

Using an instrumented hamm er, the gun barrel was tapped near its base. 
The resulting vibrations were recorded opposite the impact point with a single 
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barrel made of brick elements. 



barrel made of beam elements. 




transducer or accelerometer. Appendix A shows a plot of the recorded response for 
frequencies up to 3000 Hz. The excitation force due to the gun firing at a 
firequency of 75 Hz defines the important fiequency range for this problem. The 
normal modes of vibration of the two versions of the gun barrel were calculated in 
the free-firee configuration. Here, kinematic degrees of fireedom were used to 
simulate the restraints provided by the two low stiffiaess cords in the experimental 
set-up. Figure 5-3 shows the first mode of vibration of the advanced barrel and 
Figure 5-4 shows that of the simple barrel. Table 5-1 lists the experimentally 
determined modes of vibration of the real barrel and those calculated for the 
computer models. Figure 5-5 graphically illustrates this data. 



Figure 5-3. The first normal mode of vibration of the advanced barrel. 







TABLE 5-1. COMPARISON OF ACTUAL AND MODELED GUN BARRELS 


Modes of Vibration 




Experimental Barrel Data 

Simple Barrel Data 

Advanced Barrel Data 

72.8 

71.6 

64.8 

217 

224 

201 

431 

455 

400 

694 

744 

652 

1040 

874 

809 

1428 

1111 

982 

1848 

1515 

1357 

2363 

1550 

1501 

2913 

2013 

1775 


Looking at Figure 5-5, it can be seen that there is good correlation with both 
models and that of the real barrel over the first four modes of vibration or up 
through roughly 700 Hz. According to general finite element modeling practice, 
this agreement of the first four modes suggests that both gun models are validated. 
The agreement through 700 Hz for both g\in models is considered excellent. It is 
also interesting to note that the simple barrel data matched the experimental 
barrel data more closely at the higher frequencies than the advanced barrel. This 
was probably due to the fact that the single layer of brick elements used in the 
advanced barrel model were not representing beam behavior very well. As 







































Figure 5-4. The first normal mode of vibration of the simple barrel. 

mentioned in the previous chapter, more bricks are required to reflect bending and 
shear than do beam elements. Therefore, in the case of the advanced barrel, there 
should at least be two layers of brick elements used to define the barrel’s thickness. 

B. MODES OF VIBRATION OF THE COMPLETE GUN MODEL 

After completion of the PHALANX gun model, its static properties were 
validated against the actual gun. The next step was to calculate its normal modes 
of vibration. Table 5-2 lists the first 21 calculated modal frequencies. Upon 
observing the normal modes, it became apparent that the first mode, which had a 
frequency around 0.3 Hz, was the normal mode of vibration of rotation about the 
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Modes of Vibration 



Figure 5>5. Comparison of the modes of vibration of the real barrel versus 
the two computer models. 


lengthwise axis. This is a computational error since the gun barrel assembly is 
free to spin in the model and does not have any torsional restoring force associated 
with this motion. Therefore, it has a nonlinear, rigid body rotational mode of 
vibration of essentially 0.0 Hz. The fact that it was not exactly zero was probably 
due to the bearing springs not having coincident nodes which creates an unrealistic 
moment arm (see Chapter VII. Suggestions for Future Work). As mentioned in 
the previous chapter, the normal modes that are being calculated are those without 
gyroscopics or an unbalance (centrifugal) excitation. Therefore, the normal modes 
of the gun were resolved, this time restricting the rotational motion of the rotor, 
which removed that first torsional mode. Hence, the first 20 flexible modes were 
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calculated. The periodic force due to the firing of the gun is at a frequency of 75 
Hz and lies between modes 6 and 7. A complete pictorial record of all the modes 
used in the dynamic analysis of the gun appears in Appendix B. 

It is also important to note that in the n-^de shape pictures to follow, the 
amplitude of the deformation is greatly exaggerated in order to the general 

shape of the normal mode and the relative displacements. Correct amplitudes will 
be found later when the dynamic response of the system is evaluated after driving 
it with the periodic firing force and values of damping are assigned. Obtaining 
normal modes is an important step in the process of obtaining displacement 
information. 




TABLE 5-2. MODES OF VIBRATION OF THE PHALANX GUN 


Modes of Vibration 

Normal Mode 

Frequency (Hz) 

Normal Mode 

Frequency (Hz) 

1 

0 

12 

231 

2 

14.4 

13 

242 

3 

15.5 

14 

246 

4 

16.0 

15 

254 

5 

35.3 

16 

254 

6 

55.3 

17 

274 

7 

111 

18 

274 

8 

116 

19 

289 

9 

129 

20 

308 

10 

165 

21 

309 

11 

178 




1. First Flexible Mode 

The first mode, at a frequency of 14.4 Hz, is considered to be significant 
in the analysis of the dispersion data. This is the second normal mode in Table 5- 
2. It appears to be a link between the periodic excitation force and its rotation 
which may be causing a resulting whirling motion of the six-barrel system. Here, 
the barrels act like an overhung shaft which whips around its center line with the 
barrel tips tracing a new, larger circle. This ty]7e of motion may provide an 
explanation for dispersion patterns that appear to be a circle of the six individual 
barrel patterns as mentioned in the first chapter. This motion will be discussed in 
greater detail in the last chapter. 
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Figtire 5*6 shows a snapshot of this first flexible mode. It primarily 
consists of a rigid body motion where the barrels are straight and pitch from side 
to side. The rotor through the barrels move with respect to a stationary gun body. 
These components are pivoting about the thrust and needle bearings. There is 
little, if any, fore and aft motion. There is no recoil motion as indicated by the 
recoil adapter springs which are not exercised. In this mode and all the modes to 
follow, one component that appears to be "taking a lot of stress" is the rear gun 
body wall. Its deformation is readily apparent when a transparent picture is taken 
showing its motion with respect to the rotor, as can be seen in Figure 5-7. 



Figure 5-6. Hidden line image of the first flexible mode. 
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Figure 5*7. Hidden line image showing the deformation of the rear gun 
body plate. 

2. Second Flexible Mode 

The second flexible mode, at 15.5 Hz, is a combination of fore and aft 
motion and pitching of the gun barrel assembly. This is the third normal mode in 
Table 5-2. Figure 5-8 shows that the barrel assembly remains rigid as it pitches in 
the vertical plane. At the same time, the whole gun recoils and exercises the recoil 
adapter springs. It is also apparent that the gun body is also being exercised in 
this mode for it shows some deformation. 

This is one mode where an easy hand calculation can be done in order to 
ensure that the numerical value of the normal mode frequency is reasonable. 
Assuming that the rotor, barrel assembly and gun body are rigid, the motion can 
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Figure 5-8. Hidden line image of the second flexible mode. 


be described as a symmetric mode of the recoil adapters. That is. the whole gun 
of mass M = 120 kg is attached to ground via the two recoil adapter springs as 
shown in Figure 5-9. 

The equation of motion for such a system is 

iWi'+2K)e{Xo-x) = 0 . (5.1) 

Here, Xo is the equilibrium position of the lumped mass before it vibrates a 
distance x along the x axis. Kr = 590,485 cN/cm, is the stiilhess value for the 
recoil adapter springs used in the numerical model. Prom the equation of motion, 
the normal mode angular fi*equency of vibration, u, is 
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This compares favorably with the model value of f = 15.5 Hz. The analytic 
calculation assumed a mode of oscillation purely along the x axis, whereas in the 
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numerical model, there is a pitching motion aa well as the recoil motion which 
would contribute to the determination of its modal frequency. 

3. Third Flexible Mode 

The third flexible mode, or fourth normal mode, is similar to the previous 
mode. A snapshot of its motion is shown in Figure 5-10. Its frequency is only 
slightly greater at 16.0 Hz. The one apparent difference is that there is less recoil 
motion than seen in the previous mode. This is evidenced by the lack of 
deformation in the recoil adapter springs. Its motion is almost a purely pitching 
motion of the barrel assembly in the vertical direction. 



Figure 5-10. Hidden line image of the third flexible mode. 








4. Fourth Flexible Mode 


The first thing noticeably different in the fourth flexible mode, or fifth 
normal mode, is the massive deformation of the gun body as displayed in Figure 5- 
11. The mode of oscillation of the g\in body shows a bending node in its mid- 
section with a normal mode frequency of 35.3 Hz. The gun barrel assembly 
remains rigid with respect to the thrust and needle bearings. The motion of the 
barrel assembly appears to be caused solely by the bending of the gun body as it 
heaves up and down. There is no recoil motion in this mode. It is important again 
to note the deformation of the rear gun body plate which seems to suggest that its 
strength would play a crucial role in restricting barrel tip displacement. 



Figure 5-11. Hidden line image of the fourth flexible mode. 
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5. Fifth Flexible Mode 


The fifth flexible mode, or fourth normal mode is at 55.3 Hz. Again it 
appears to involve a mode of oscillation of the gim body alone as shown in Figvire 
5-12. However, the gun body oscillation is that of a breathing mode where the 
walls move in and out. Like the previous mode the gun barrel assembly remains 
stationary with respect to the bearings. So the barrel motion is directly related to 
the motion of the gun body walls. As the walls breath, the recoil adapters appear 
to twist in the vertical direction. There is no motion in the fore and aft direction 
except that caused by gun body oscillation. 



Figure 5-12. Hidden line image of the fifth flexible mode. 
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6. Sixth Flexible Mode 


The final mode to be discussed is the seventh normal mode and is the 
first mode in which the barrels themselves are exercised. This mode, at 111 Hz, is 
shown in Figure 5*13. The barrel tips remain parallel to ground as they oscillate 
about a single bending node located just forward of the mid-barrel clamp. The 
boundary conditions of the muzzle clamp seems to be vital in relation to barrel tip 
motion. In the real g^, three of the barrels are unlocked in the muzzle clamp 
while the other three are locked. In this model, all the barrels are rigidly attached 
or locked to the muzzle clamps. A possible future modification would be to unlock 
three of the barrels in the muzzle clamp and simulate the effect of a gap around 
them using spring elements. The difficulty would come in determining the stiffness 
values of the springs since this is a nonlinear problem. A new version of I-DEAS, 
called the Master Series, is able to simulate the effect of a gap or tolerance. After 
the real g\in is set up on-sight, then tests can be performed to determine the 
proper stiffness values. This vibrational mode of the barrels with clamps seems to 
exhibit the behavior of a deep flange beam which is different than the six barrels 
represented as a single beam. A deep flange beam tends to shear like the barrels 
in the model rather than bend like a beam. 

Looking at the rest of the gun, there is massive deformation of the gun 
body with a single bending mode in the middle. The rotor essentially remains 
parallel to ground as the gun body and barrels flex. There is no recoil motion 
associated with this normal mode. 

C. MODAL TRUNCATION 

With 22,470 degrees of freedom, there are an equal number of normal modes 
of vibration. Solving for all of these normal modes is unnecessary, but would 
defeat the purpose of reducing the complex model into relatively simple normal 
modes. The System Dynamics Analysis module, to be discussed fiirther in the next 


86 





Figure 5*13. Hidden line image of the sixth flexible mode. 

chapter, calculates the dynamic response of a system by first solving for the modes 
of vibration. Instead of using the 22,470 physical degrees of freedom, it can 
describe the response using the first 20 flexible modes of Table 5-2. 

According to generally accepted finite element modeling practice, the number 
of normal modes that should be carried through for a dynamic response analysis 
are those that fall within twice the excitation firequency. Since the gun’s principal 
excitation fiiequency is 75 Hz, the modes through 150 Hz should be calculated. By 
solving for the first 20 modes, up to a frequency of 309 Hz then this criteria i 
more than satisfied. 
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VI. DRIVING AND DAMPING OF THE PHALANX GUN 


A. DYNAMIC ANALYSIS OF THE GUN MODEL 

Having completed a modal analysis of the PHALANX gun model, the next 
step, as outlined in Chapter III, is to perform a forced response analysis of the 
gun. The excitation forces caused by the bullets recoil and estimates of damping in 
the gun are applied. The resulting response, such as barrel tip displacement, can 
then be determined and compared to experimental data. The most prominent 
normal mode can be found to describe the actual motion of the gun when firing. 
With the completion of these steps, a working model of the PHALANX gun has 
been created. The factors which contribute to dispersion can be studied 
systematically, and potential solutions to problems can be tested to provide an 
additional source of information to the PHALANX program office. Various 
components can be "swapped" out to evaluate modified components, and the 
resulting effects on gun performance predicted prior to any hardware tests. This is 
the first finite element model of its kind ever constructed of the PHALANX gun, 
and the possibilities for its use are only limited by imagination. 

The I-DEAS computer code evaluates a structxire’s response to applied forces 
in the computer module called System Dynamics Analysis (SYSTAN). The method 
used by SYSTAN is to first develop a modal component model of the structure. 
That is, a certain number of normal modes of a structure are calculated and stored 
as a new modal model of the structure which is then used to calculate the dynamic 
response. Other methods can be used, such as directly integrating the equation of 
motion 


F = 


M 


A' + 


X + 


K 


( 6 . 1 ) 


Here, column matrices are represented with underlines and square matrices with 
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brackets. F represents the applied forces and their direction on each degree of 
freedom in the system. [M] is the modal mass matrix that was obtained through 
the process of transforming to normal coordinates described in Chapter III. Each 
of the ith matrix elements is a measure of how much mass is participating in the 
ith normal mode. Similarly, [K] is the modal stiffness matrix and is a measure of 
the stiffness occurring in each of the normal modes. [C] is the associated damping 
for and between each degree of freedom. The displacements, velocities and 
accelerations for each degree of freedom are represented by X ,X , and X 
respectively. However, when there is a structure with many degrees of freedom, 
like the PHALANX gun model with 22,470, it becomes apparent that the 
calculations required are memory and time consuming. Reducing the number of 
physical degrees of freedom, X, &x>m 22,470 to 20 modal degrees of freedom, x, 
greatly reduces the complexity of the model. In order to accomplish this reduction, 
a transformation to normal coordinates diagonalizes the mass, stiffness and 
damping matrices of the equation of motion as discussed in Chapter III. 
Reductions which do not diagonalize are also possible, but will not be discussed 
here. The equation of motion can then be written as 


Mu 0 . . ■ 


0 


Kn 0 . 

0 M 22 0 


0 C 22 0 


0 K 22 0 

0 . 0 

X + 

0 

0 

X + 

0 

0 

. 0 My 


• 0 Cy 


. 0 Ky 


This equation is similar to the diagonalized equation of motion (3.60), developed in 
Chapter III, except that now a damping matrix [C ] and an excitation force F have 
been added. Damping is necessary when evaluating a dynamic response in order 
to simulate the effects of energy dissipation in each normal mode. Otherwise, the 
motion of a structure will continue to grow unbounded as the excitation force is 
applied. In (6.2) and in SYSTAN, the damping is assumed to be viscous damping 
proportional to the velocity of the structure simulating the effects of friction. With 
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the matrices diagonalized, and hence uncoupled, the equations of motion are easily 
solved. [Ref. 11] 


The conversion can then be made back to the physical deg^rees of freedom 
through the relation 


^ = 



(6.3) 


where [< 1 >] represents the spatial and temporal parts of the mode shape as 
discussed in Chapter III. Each column of this matrix is an eigenvector 0 ‘ which 
corresponds to each modal degree of freedom. [Ref. 11] 


B. DAMPING IN THE PHALANX GUN 

Damping in the PHALANX gun model is defined as Rayleigh or proportional 
damping. Since the diagonalized modal mass and stiilhess matrices were obtained 
without the velocity damping term iu the equation of motion, then the damping 
matrix [C] is not diagonahzed. Instead, the damping coefficients are defined by 
the relation 

CsoJf + pAf (6.4) 

Diagonalizing this matrix yields 

0 ^C0 = + p0^Af0 (6.5) 

= a[K] + p[Ml , 

where [K] and [M] are now diagonalized matrices. Substituting these relations 
back into the equation of motion gives 

F = [M]i + [K]x + (a[K] + p[M])x . ( 6 . 6 ) 

The quantity in parenthesis is called modal damping and it was this quantity that 
was assigned a nximerical value in the gun model. The correct modal damping 
values, a and P, can only be determined by measurements on the actual system 
being modeled. Therefore, damping is associated with each normal mode and not 
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with each element in the gun model. To more accurately represent damping of the 
real gun, damping would have to be defined in terms of two physical points, such 
as a dashpot in the recoil adapter, and not just with each mode. Each mode itself 
is not a pure motion as well. For example, the recoil mode has a certain amount of 
pitching associated with it, so an additional value of damping should be assigned to 
this additional motion. In the absence of measurements, the particular 
components of the gun involved in each normal mode can be iised to estimate 
dam ping. Damping values were assigned based on knowledge of how the 
components are assembled, for example, welded or bolted components have certain 
damping, while the barrels alone have only small damping as observed in the 
modal vibration tests and the recoil adapter springs are intentionally damped by 
the manufacturer, as described by a hysteresis plot. 

In the SYSTAN module, damping is defined as a percent of critical damping. 
To clarify further, the effects of damping are associated with three physically 
distinct cases. These cases are derived firom the equation of motion (6.2) where an 
exponential function of the form 

X=e’‘ , (6.7) 


was used as a solution with the driving force set to zero. The equation of motion 
for the first mode is then 


+ Kue’* = 0 . (6.8) 

dt^ dt 

Performing the necessary derivatives and canceling the factors of yields the 
quadratic equation 

+ Ciig + Kii = 0 (6.9) 


Solving for q using the quadratic formula gives 

X 

_ -C|i ± (Cii^ - 4M^iKn)^ 

■ 2Mh 


( 6 . 10 ) 


91 





The three difTerent damping cases are then determined from the value of the 
discriminant Cn^ - 4MiiKii. Ouerdamping occurs for the case where 
> 0. In this case there is an exponential type of decay. Given an 
initial displacement, the structure does not oscillate about its equilibrium position 
as shown in Figiire 6-Kb). The critical damped case occurs when 

Cii^ - = 0. Like overdamping, critical damping does not oscillate, but 

dampens out asymptotically as shown in Figxire 6-l(c). Finally, the underdamped 
case occurs when Cn^ - 4MiiKii < 0 where the structure oscillates about its 
equilibrium point sinusoidally, but eventually dampens out as shown in Figirre 6- 
1(a). [Ref. 18] 


THE DAMPED HARMONIC OSOULATOR 


t 



Remm of hamioaic oKillitor to eqnUibiinin. (a) Underdampcd. (b) Ovudampod. 
|c) CritiaU; damped. 


Figure 6-1. Graphs of displacement for the three case of damping. 
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When determining a damping value based on knowledge of how the 
components are connected, experts in the industry have accepted general empirical 
values until more information is obtained. Modes involving welded parts are 
assumed to have damping values of around 1 percent of critical damping whereas 
bolted structures have values between 3 and 5 percent. For clamped components, 
2 percent is an accepted value. Since most of the components are clamped in the 
gun, 2 percent was used for most of the modes. [Ref. 11] 

A second case estimates the modal damping value for a single gun barrel. 
Only this component has been observed while measuring the modal vibration. 
While obtaining the normal modes of vibration of the single g\jn barrel by tapping 
it with the instrumented hammer, it was observed to have a very low damping 
value. It was observed to have a damping value of less than 1 percent in the free- 
free configuration. The damping value used in the computer model was raised to 2 
percent to reflect the fact that each gun barrel is clamped at its base, middle and 
end points. 

Another component that has experimental damping information available are 
the two recoil adapters. A hysteresis plot provided by Martin Marietta for the two 
recoil adapters is shown in Figure 6-2 [Ref. 15]. It shows the loading force applied 
to the recoil springs for one cycle on the vertical axis and the resulting deflection on 
the horizontal axis. The upper branch reflects the initial loading and the lower 
branch the return to battery. The area contained within the two branches is the 
energy dissipated per cycle, £/)= 4.34 x 10^N-cm (384 lb-in). The damping 
coefficient, Cfty,t, is then determined from the equation 

Ed = , ( 6 . 11 ) 

where, x is the amplitude of oscillation about its driven equilibrium point 
estimated by the mid-point in the hysteresis plot, or x = 0.4 cm (0.15 in). After the 
first few rounds are fired, the gun sets back about 1.3 cm (0.5 in), and then 
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oscillates about this new equilibrium position xintU the gun stops firing. It never 
returns to battery between shots Like other naval guns. The angular frequency of 
oscillation of the recoil mode O) = 30it/ s. Rearranging and solving for yields 


Chytt 


—= 91.6 
jj be ro) cm / s 


( 6 . 12 ) 


The dam ping coefficient, y= Ctty^/ Ccru> used by the computer code is expressed as 
a percentage of critical damping, Cent where 

Celit = AmK , ( 6 . 13 ) 


m is the mass of the gun, and K is the stiffness of the recoil adapters. 
Substituting values into this equation gives 


Ceru = 2.38 X 10^- 


N 


cm / s 


( 6 . 14 ) 


so that y can now be calculated as 


Y=.^ = 0.39 

^erit 


(6.15) 


The value for y in (6.15) is consistent with the actual damping value used in 
the design of the recoil adapters, y = 0.33, as provided by Martin Marietta [Ref. 
19]. A damping coefficient which approaches critical damping, y 1, gives a more 
optimum return to equilibrium after the excitation force has stopped. In the 
computer model a damping value of y= 0.33 was used for this mode. The recoil 
mode, referred to in Chapter V as the third flexible mode of vibration at 15.5 Hz, 
shows the gun sliding fore and aft exercising the recoil adapter springs with a 
slight pitching of the barrels in the vertical plane. This is the only mode using a 
damping value other than 2 percent. The purpose of the recoil adapters is to 
dissipate some of the force caused by the firing rounds and to couple the remaining 
force to the ship. Allowing a gun to recoil even a few calibers greatly reduces the 
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Figure 6-2. Hysteresis plot of the recoil adapter springs. 

force of recoil transmitted to the mount so that a much Ughter gun mount can be 
used. 

Table 6-1 lists the damping values applied to each normal mode of vibration 
used in the PHALANX gun model. Most of the damping values used in the model 
are crude estimates based on information derived without the actual gun on sight. 
Further experimental tests should be performed allowing for further refinement of 
the damping values. 


95 




TABLE 6-1. MODAL DAMPING VALUES 


Damping of the Modes of Vibration 

Normal Mode 

Damping (y) 

Normal Mode 

Damping (y) 

1 

0 

12 

0.02 

2 (yaw) 

0.02 

13 

0.02 

j 

3 (recoil, pitch) 

0.33 

14 

0.02 

4 (pitch) 

0.02 

15 

0.02 

5 (gun body mode) 

0.02 

16 

0.02 

6 (gun body mode) 

0.02 

17 

0.02 

7 (barrel bending) 

0.02 

18 

0.02 

8 

0.02 

19 

0.02 

9 

0.02 

20 

0.02 

10 

0.02 

21 

0.02 

11 

0.02 




C. DRIVING THE PHALANX GUN 

The principal excitation force of the PHALANX gun is the repeated firing of 
the rounds. Only one barrel shoots at a time, so that a firing rate of 4500 shots 
per minute gives the frequency of the excitation force 75 Hz. The experimentally 
determined force on the breech for a single shot lasts roughly 2.3 milliseconds as 
displayed in Table 6-2 and plotted in Figure 6-3 [Ref. 15]. 
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TABLE 6-2. FORCE ON GUN BREECH 


Force of Single Shot on Breech 

Time (milhseconds) 

Force (lb) 

Time (milhsecond) 

Force (lb) 

0.05 

8.74 X 10^ 

0.80 

2.85 X 10“ 

0.10 

1.26 X 10^ 

0.85 

3.02 X 10“ 

0.15 

1.72 X 10^ 

0.90 

3.12 X 10“ 

0.20 

2.33 X 10^ 

0.95 

3.16 X 10“ 

0.25 

3.12 X 10^ 

1.00 

3.13 X 10“ 

0.30 

4.14 X 10^ 

1.05 

3.02 X 10“ 

0.35 

5.45 X 10^ 

1.15 

2.59 X 10“ 

0.40 

7.11 X 10^ 

1.25 

2.02 X 10“ 

0.45 

9.17 X 10^ 

1.45 

1.16 X 10“ 

0.50 

1.17 X 10“ 

1.65 

7.68 X 10^ 

0.55 

1.45 X 10“ 

1.85 

5.57 X 10^ 

0.60 

1.76 X 10“ 

2.05 

4.29 X 10^ 

0.65 

2.08 X 10“ 

2.25 

3.46 X 10^ 

0.70 

2.37 X 10“ 

2.31 

3.26 X 10^ 

0.75 

2.64 X 10“ 




A t 3 rpical missile engagement burst lasts for only a few seconds up to a 
maximum of 5 seconds before the barrels deteriorate due to the high temperature. 
Therefore, a driving force of one hundred shots was constructed. Due to the 
sharpness of the firing "spike" compared to the length of time between shots, 13.3 
milliseconds at 75 Hz, it was necessary to determine more data points for each 
shot in order to increase the resolution of the response. A function was fit to the 
available data in Figure 6-3, as shown in Figure 6-4. The function’s representation 
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Figure 6-3. Plot of the force on the breech from a single shot [Ref. 15). 
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of the first 3 shots is shown in Figure 6-5 with dots representing the experimental 
data points described above. The function was used in SYSTAN to create a forcing 
function using 26,400 data points over a one hundred shot burst. The 26,400 
evenly spaced points are spaced at intervals of 0.05 milliseconds to provide a good 
representation of each 2.3 millisecond shot. The number of data points used is 
then determined as 


Number of points 


- ^-3 seconds _ 26,400 

0.05 milliseconds 


where, 1.3 seconds is the length of time for a 100 shot burst. 


( 6 . 16 ) 
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Force on Breech 
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Figure 6-4. Plot of the function used to approximate the single shot data. 

A single shot hy itself is broad band and is comprised of many frequencies. It 
is Uke an impulse or delta fimction. The firing of a single shot by the gun model, 
would be expected to excite many normal modes of vibration over a wide frequency 
range. However, the gun is never used in the single shot mode, so this case is 
irrelevant. A long series of shots is no longer broad band, but has a fimdamental 
frequency of 75 Hz. This type of excitation force would be expected to excite fewer 
modes in a small range around 75 Hz. 

A burst of shots is applied at the base of the firing barrel located 
approximately between 10 and 12 degrees from top dead center. The firing 
sequence was first designed for the older firing rate of 3000 shots per minute. In 
this case, the round was ignited at this off-center position so that when the barrel 
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Figure 6-5. The function’s approximation for the hrst 3 shots. 


reached the top, the round would just be leaving the barrel. However, the firing 
rate has been changed to 4500 shots per minute without an adjustment being 
made to the igniting point. Because of the higher rotation rate of the barrel 
assembly, the bullet leaves the barrel approximately 14 degrees past top dead 
center. The time varying magnitude of the firing force is applied to the gun over a 
20 degree arc. To simulate this off-center firing force, the barrels were moved 10 
degrees from the top. The gun model barrels can be rotated to any point along this 
20 degree arc, and the fining force reapplied. The barrel tip motion can then be 
recalculated. This is an approach being taken by Martin Marietta. It has not been 
possible to continuously rotate the barrels in the numerical model while applying 
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the firing force. The current arrangement gives a good approximation. It appears 
that there is a way to "trick” the computer code into simulating the effects of 
rotation, but this is ctirrently beyond the scope of the airrent work. Exploring this 
improvement is recommended as a further refinement of the model. 

A possible comphcation in applying the firing force to the gun is the fact that 
it may be a follower force problem. That is, as the barrel assembly becomes 
displaced from its center line, for whatever reason, the firing force ought to be 
applied to the gun in this displaced position. The combination of this displacement 
and normal firing would add other, currently unpredictable, effects on the barrel 
tip motion. The following force problem is probably best suited to be described 
analytically. Preliminary steps in this direction have been taken by Professor 
Gordis. 

Another force acting on the gun arises firom the rotating barrels. The rotation 
rate of the six-barrel system is 12.5 Hz = 75 Hz/ 6 barrels. This motion can play 
an important role in gtm dynamics since there is a nearly resonant normal mode of 
vibration at 14.4 Hz. This normal mode describes a pitching motion of the barrels 
and may contribute to the dispersion pattern if excited. 

When the rotating barrels are displaced firom their center line, they will cause 
excitation. A perturbing force that could cause a displacement would be a weight 
imbalance between the barrels. The mass of a barrel may vary by as much as 
0.227 kg (0.5 lb) away from the average weight of 8 kg (18 lb). This imbalance 
would cause the barrel cluster to bend slightly as they rotate. This displacement 
could also start when the barrel cluster kicks up as the gun first fired. The kicking 
up is caused by the placement of the recoil adapters and the rear ball joint in 
relation to the point of application of the firing force. Displacement is created due 
to the imbalance of the firing point and the three supports. Any of these causes 
leads to excitation of a physical motion called gyroscopic "whirl." The near 
resonance of the barrel rotation rate and the first flexible mode of vibration 
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appears sufficient to excite this pitching mode and greatly edTect dispersion. 
Furthermore, a non-preloaded thrust bearing, as is used in the PHALANX gun, is 
considered to be a root cause of gyroscopic whirl [Ref. 20]. Gyroscopic whirl will be 
further discussed in a latter section. There are a variety of "causes ' for dispersion 
that may feed upon others. If one could determine the primary causes then it may 
be possible to prevent or limit others. 

Another excitation force on the gun is the torque exerted on the barrel as 
each round spins up as it travels along the barrel. The barrel lands have a right- 
hand exponential gain twist which rotates the penetrator to a rotation rate of 
120,000 revolutions per minute (2,000 revolutions/ second) as it leaves the barrel 
[Ref. 21]. The torque to achieve this rotation rate can be calculated by first 
determining its rate of angular acceleration. 


t. 


(6.17) 


where, (O/ is the final rotation rate, % is the initial rotation rate, and is the 
time that the round travels the length of the barrel. Substituting the values in 
yields 

i, = 5.46 X 10'. 

second^ 

The torque applied by the barrel on the round to achieve this rotation rate can be 
calculated using the following relation. 



(6.18) 


where I is the moment of inertia of a cylinder which approximates that of the 
roxind, and a is the inner radius of the barrel. I is determined according to the 
relation 


frirO^ 

2 


(6.19) 
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where is the mass of the round. / is calculated to be 

, ^ (0.071 tgxio mm-f ^ 3 x 10-» , 

A 

Substituting this value into equation (6.18) gives 

Fft = 435.6 lb , 

showing a significant force on the barrels. This type of force may be difficult to 
simulate with this model due to its variation in location with time, and is similar to 
the rotating barrel problem. 

Other possible forces include an imbalanced load on the gun with heavy 
preshot rounds on one side and the li|^ter expended casings on the other. Also, 
the act of loading a round is not insi gnifi cant and would add an impulsing torque 
to the gun assembly. Furthermore, the actual rate of firing varies since the 
hydraulic motor which rotates the harrel assembly, and determines the firing rate, 
varies its speed by as much as ± 50 revolutions per minute. These variations could 
occur as sudden rotation speed shifts resulting in a torque or "jerk" on the barrel 
assembly. 

The modeling of the PHALANX gun in this study ends at the three 
supporting fixtures. Other factors affecting the gun’s motion from beyond these 
points are not being simulated, but may also play a role in the actual dispersion. 
The gun mount will rarely be in a stationary position when attempting to shoot 
down targets moving at speeds on the order of a few Mach. The gun can slew at a 
maximum rate of 92 degrees per second in elevation and 126 degrees per second in 
train or rotation. This motion, may serve as a perturbing force to start the 
gyroscopic whirl mentioned earlier. The whole gun mount is mounted on a set of 
shock absorbers whose purpose is to limit the shock wave of a close-order 
explosion. The vibrating gun interacting with these shock mounts may merit 
future consideration. 
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Ib an actual gxin, all of these forces may play a role in the observed dispersion 
pattern. Some of these forces may be added to the current gun model analysis, but 
this study only addresses the role played by the most significant force of the gun 
firing. 

D. RESPONSE OF THE GUN MODEL 

While driving the gun model with a one hundred shot forcing function, the 
resulting response of the top barrel tip was evaluated to its translation along each 
of the three coordinate axes. Shown in Figure 6-6 is the amplitude of the barrel 
tip’s motion along the x-axis in the fore-aft direction versus time. The plot shows 
that the gun immediately sets back by roughly 1.27 cm (0.5 in) after the first shot 
is fired. It then begins to oscillate about its new eqiiilibrium position by roughly 
0.80 cm (0.32 in) at a firequency between 15 to 20 Hz as one or more of the lower 
order modes is excited. It appears to reach a steady-state oscillation near the end 
of the burst of shots. When the firing is complete after 1.3 seconds time on the 
horizontal axis, the gun moves forward to its initial position and vibrates like an 
underdamp)ed harmonic oscillator. During excitation, the displacement versus time 
shows an irregular pattern of jagged peaks indicating excitation of several normal 
modes. Understanding which normal modes are involved in the response will 
provide the greatest benefit in predicting the dispersion pattern. The 
determination of these modes will be provided by a frequency analysis which will 
be the subject of the next section. 

The important barrel tip displacement in the vertical direction is shown in 
Figure 6-7. The barrel tip can be seen to pitch upward significantly after the first 
shot, and then approach steady-state oscillation near the end of the firing burst. 
The peak-to-trough barrel tip displacement measured near the steady-state motion 
is about 2.5 mm (0.1 in) and gives the observed 1.64 milliradian dispersion for a 
1.52 m long barrel. The first 20 to 30 shots of the firing burst have a much greater 
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Figure 6-6. Response plot of the barrel tip in the fore-aft direction 
(x-plane). 

dispersion pattern than the last 20 to 30 shots. No data is currently available for 
comparison to this modeled phenomenon. Since the recoil normal mode of 
vibration at 15.5 Hz has an associated pitching motion, it appears that this motion 
is important for predicting the dispersion pattern. Unlike the response in the 
fore-aft direction, this vertical response appears to involve only a single mode with 
a h quency of oscillation between 15 to 20 Hz. The barrel tip displacement 
corresponds to an angle of 1.64 milhradian movement which represents a 
significant start in quantifying the amount of dispersion associated with barrel tip 
motion. 
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Figure 6-7. Response plot of the barrel tip in the vertical direction 
(y-plane). 

Figure 6-8 shows barrel tip translation along the z-axis, or side-to-side 
motion. The amplitude of motion is noticeably less than in the vertical plane and 
demonstrates the relative strengths of the moment arms about the center of 
gravity. In the vertical plane, there is a larger moment arm from the firing point 
to the recoil adapter pivot points than there is in the z-plane so that greater 
displacement in the y-p’-ne is expected. The off-center firing force causes a larger 
displacement in one direction than in the other. This shows that actual dispersion 
should be broken down into components in the two transverse dimensions x and z. 
Based on this analysis, the Naval Surface Warfare Center, Dahlgren, and Martin 
Marietta have now begun to describe their observed dispersion patterns as 
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elliptically shaped [Ref. 22]. The maximum barrel tip displacement in the z 
direction is about 5.10 mm (0.20 in), and settles into a steady-state displacement 
of about 2 mm (0.08 in). The response in the z-plane appears to be multimode 
which is similar to that of the x-plane response. The response and ring down show 
beats which suggests that the multimode response involves two normal modes 
whose frequencies are close together. Recall that the response of a structure is just 
the summation of the excited normal modes of vibration as given by 


X = 



( 6 . 20 ) 


The beating appearing in this response plot indicates that the normal modes 
excited have similar damping. The multiple normal modes excited in the x-plane 
apparently had different damping. More will be discovered in the next section on 
this topic using the firequency response plots. 


E. FREQUENCY RESPONSE FUNCTIONS 

The previously discussed response plots are important for ensuring that the 
hnite element model appears to be a close approximation of the real structure. It 
is a check on the assvimptions made during the modeling process. However, it is 
more important to understand the normal modes involved in the response. In this 
way, it is possible to pinpoint causes of large amplitude displacements of the 
structure. A structure’s response to an excitation depends on the magnitude of the 
excitation force, the frequency of the excitation force, the point of application and 
damping. Large displacements can occur if the firequency of the excitation is near 
that of a normal mode of vibration. This process is called resonance. Here it is 
possible to get large structural displacements even though the magnitude of the 
excitation may be small. If the magnitude of the driving force is large, then the 
structural deformation will also be large. It is also important whether the point of 
application of the driving force is at a node of the mode shape or at an antinode. If 
applied at a node the response is less than if the force was applied at an emtinode. 
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Figure 6-8. Respionse plot of the barrel tip in the horizontal direction 
(z-plane). 

Therefore, finequency response functions and modal analysis are crucial in 
diagnosing the barrel tip displacement of the PHALANX gim. 

One of the best tools for understanding which normal modes are involved in a 
forced response is called the frequency response function (FRF). This function is 
obtained by applying an input force at one point in the structure and measxiring 
the response such as displacement at another point. The FRF is then the ratio 
the response divided by the input. In this study, three FRFs were obtained, one 
for response of the barrel tip motion along each coordinate axis. Mathematically, 
they appear as 

•^3748^0}) Y 3748(Cfl) (6 21) 

F, (w) ’ F, (w) ’ F, (o) ^ ' 
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In these three functions, the numerator stands for the displacement in the 
respective direction at node number 3748 which is at the tip of the bring barrel. 
The denominator represents the driving force function along the x (fore-afl) axis 
applied at node number 1165 at the base of the firing barrel. To determine the 
FRF, the firequency of the force is swept from zero to a high value. The FRF 
appears as a series of spikes which are resonant with the normal modes of 
vibration. The height of the spikes depends on the magnitude of excitation and the 
width of the peaks depends on the damping associated with that normal mode. A 
wider spike has a larger damping value. This FRF is complex since both 
amplitude and phase of oscillation are involved. [Ref. 11] 

Figure 6-9 is a FRF of the gun response along the x direction at the barrel 
tip. As expected finom the transient plot in Figure 6-6, there is a multimode 
response involving the first flexible mode at 14.4 Hz and second at 15.5 Hz which 
are being excited. The second mode is much more excited than the first based on 
the height of the two peaks. This is expected since the second mode is the recoil 
mode, while the first is a rigid body pitching mode. It appears also that no other 
modes are being excited. The next three or four modes involve the deformation of 
the gun body, and do not appear to be excited. This may be because the point of 
appbcation of the driving force is at the rotor which is suspended by the recoil 
bearings. 

The next FRF is shown in Figure 6-10 describing the barrel tip response in 
the y direction. It appears to verify the transient response in that this is a single 
mode excitation. The mode that appears to be excited is the recoil mode mixed 
with some pitching motion. This is identified as the second mode due to the fact 
that its normal mode motion is in the vertical plane, and that the spike is fairly 
wide suggesting that it uses the higher damping value of 0.33 rather than the 
others at 0.2. 
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Figure 6-9. Frequency response function of the barrel tip in the x direction. 

The last FRF shown in Figiire 6-11 is that of the barrel tip in the z direction. 
As mentioned earlier, the response plot of Figure 6-8 showed that the phenomenon 
of beating was occurring. Studying the FRF, it appears that the first flexible mode 
at 14.4 Hz and third flexible mode at 16.0 Hz are being excited. These two modes 
have similar damping values allowing them to interfere and beat. Beating occurs 
when two modes are closely spaced in firequency and amplitude as seen here. The 
amplitudes were different in the x direction case, so that beating could be seen. 
The two modes involved in the z direction FRF are more closely excited to the same 
amplitude than in the x direction where only one mode was predominantly 
excited. 
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Figure 6-10. Frequency Response Function of the barrel tip in the y 
direction. 

F. OTHER FACTORS AFFECTING DISPERSION 

In addition to the resonance described above, other physical phenomenon 
such as gyroscopic whirl and rotor imbalance may contribute to dispersion. 
Gyroscopic whirl is a self-excited vibration that occurs when there is a tangential 
force normal to the axial direction of a rotating shaft. At some critical speed, the 
whirl instability will overcome any stabilizing forces, such as those provided by 
bearings. The amphtude of the whirl motion will only be limited by nonlinear 
stabihzing forces. This type of motion is common in machinery that has an 
overhanging shaft such as that shown in Figure 6-12. In this figure, the two rear 
supports can be thought of as the rear needle bearing and the forward thrust 
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Figure 6-11. Frequency Response Function of the barrel tip in the z 
direction. 

bearing. The mass of the disc at the end of the shall represents the gun barrels 
and clamps. The PHALANX gun system is similar to that of an overhung shall. 
This fact combined with the knowledge that the thrust bearing is not preloaded are 
two key requirements for the onset of whirling. The principal excitation force of 
the firing gun induces several moment arms which would lead to forces normal to 
the axis of rotation. These indications coupled with the data on PHALANX 
dispersion, such as the pattern shown in Figure 6-13, strongly implies that 
whirling is occurring. The speed of the whirl is usuedly close to that of the rotating 
shall, and accounts for the fact that each barrel appears to shoot in the same 
direction on each cycle. Discussions with Stan Wilson of Port Hueneme indicates 
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that field experience supports this possibihty [Ref. 23]. When an old thrust 
bearing was replaced with a new one, the dispersion pattern improved 
substantially. This would seem to indicate that the stabilizing forces of a new 
thrust bearing can limit the onset of whirl [Ref. 20] 



Figure 6-12. Representation of a hung shaft demonstrating the concept 
of gyroscopic whirl 

Rotor imbalance is another phenomenon found in rotating machinery which 
could lead to barrel tip displacement. If the principal axis of rotation does not 
coincide with a fixed line in space, then wobbling may occur. This rotational 
instability can be caused by machining tolerances or mass asymmetries. Once 
started it can aggravate other sources of vibration when rotational speeds are high, 
as in the PHALANX gun. The unbalance is significant since it determines the 
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Figure 6-13. A dispersion pattern of the PHALANX gun. This particular 
pattern leads to assuming that whirl is occurring. 

degree to which any bending or flexural mode of vibration is excited. At a critical 
speed, corresponding to a resonant frequency of the system, the rotor modal 
response becomes a maximum. When a rotor spins at a speed near its critical 
speed, it tends to adopt the shape of the normal mode that is excited. The amount 
of rotor displacement depends on damping in the system as well as the modal 
components of the imbalance. This modal component of imbalance may be reduced 
by adding a nximber of correctional tuning masses that would shift the normal 
mode frequency away from a critical speed. 
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Vn. SUGGESTIONS FOR FUTURE WORK 

A. REFINING THE GUN MODEL 

The current PHALANX gun model is a first cut at attempting to understand 
the gun’s complex dynamic behavior. It was decided early that it was better to go 
through the complete simulation cycle from finite element modeling to calculating 
its dynamic response, rather than getting caught up in trying to make a perfect 
model without ever calculating its response. In this way, the whole process was 
validated and demonstrated to be a useful tool in iinderstanding the dispersion as 
well as modifying various components to see how the real gun could be improved. 
Now that the modeling process is understood finom start to finish, and valuable 
experience obtained along the way, certain recommendations stand out as the next 
logical steps in continuing this research. 

The first recommendation is to refine the gun model. The purpose of this 
refinement is primarily to reduce the number of degrees of freedom so that the 
model would have a smoother animation picture and a quicker solve time. The 
model should be simplified while retaining pertinent characteristics. This would 
serve to better illustrate the dynsmiic behavior of the gun’s mode shapes and lead 
to more insight as to how the dispersion pattern is being formed. Currently, the 
animation picture is slow and lurches from one frame to the next. The second 
purpose behind further refinements was to choose better elements to capture even 
finer details of the gxin’s behavior. Now that the real gun is better understood, 
some components of the gun should have elements that model bending and 
torsional behavior with beams and thin shells whereas others should have bricks to 
model elasticity without any preferred direction. Neither of these elements can 
perform both functions. Next, the gun’s boundary conditions are now better 
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understood; some barrels in the muzzle clamp are constrained while others are 
unconstrained. Afler attending the barrel optimization conference in Burlington, 
Vermont, May 26th and 27th, 1993, the consensus was to firmly clamp all the 
barrels. This modification could easily be evaluated with the current gun model. 
Finally, some of the currently modeled components, such as the recoil adapters, 
could be more accurately represented. 

1. Muzzle and Mid-Barrel Clamps 

In an effort to reduce the degrees of fireedom of the model, the muzzle 
clamp could be further simplified. Currently the muzzle clamp shown in Figure 3-9 
has brick elements representing the plates and thin shells representing the webbed 
meshing connecting the two. This arrangement is good for representing the 
inertial properties of the clamp, but the bricks themselves do not add much to the 
model other than increasing its complexity. Since it is unlikely that clamp shear 
plays a significant role in the gun’s dynamic behavior, such elements that model 
shearing well are not necessary. Bricks, beams and shells would be important if 
this model was used to understand the gun’s thermodynamic behavior, but that is 
c\irrently not the case. Instead it is proposed that the bricks be replaced by a 
lumped mass located at the center of gravity of the muzzle clamp plates. This 
lumped maiss is an "element" choice in I-DEAS where a structure’s mass and 
inertia can be defined at a geometric point. The lumped mass can then be 
connected to the barrels via rigid elements. Figure 7-1 shows this proposed 
arrangement. These rigid elements are massless and do not flex, but they are 
effective in attaching two structures together. In this case, the barrels are 
attached to the liimped mass. This will serve the purpose of defining the 
appropriate boundary conditions. Some barrels could be connected while others 
are not in order to model the current gun in the fleet. Alternately they all could be 
easily connected together to test the proposed gun design. What is lost in this 
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representation is the bending behavior of the plates. It is uncertain if this 
behavior is important to model or not, but experiments could be performed on the 
real gun at the Naval Postgraduate School to resolve this question. In a separate 
finite element model of the muzzle clamp, the first flexible normad mode of 
vibration occurred at about 500 Hz and is far fi-om the excitation force frequency of 
75 Hz. It is unlikely that this mode is being excited. If this did not turn out to be 
the case and the mode is excited, thin shells could be used to represent the muzzle 
clamp plates vice lumped masses and rigid elements. Few degrees of freedom 
would be added to the model, since they would still be "cheaper" to use than bricks 
and have the benefit of bending behavior. The mid-barrel clamp would be modeled 
like the muzzle clamp plate. 

What is important about the muzzle and mid-barrel clamps ^s how they 
determine the boxmdary conditions on the barrels. Currently three of the barrels 
have their rotations taken out at the muzzle clamp while the other three are 
allowed to rotate. At the mid-barrel clamp the barrels are rigidly attached. It has 
been proposed that the PHALANX gun would be improved if its barrels were all 
rigidly constrained at the muzzle clamp like that of the 30 mm GOALKEEPER g\in 
[Ref. 21]. This particular gun has a smaller dispersion pattern than the 
PHALANX gun. Therefore, it is not so important that the clamps are geometrically 
correct, but that they are d)nnamically correct. 

2. Gun Body Housing 

The next area of finite element refining is that of the gun body housing. 
Its current representation can be seen in Figure 3-€. The model could be improved 
if there was better matching of the degrees of freedom between the thin shells of 
the gvm body and the bricks of the thrust bearing. The thin shells have six degrees 
of fireedom while the bricks have three. This creates a possible boundary condition 
problem and appears to manifest itself in the extreme distortion of the gun body in 
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Figure 7-1. Recommended modification for representing the muzzle 
I clamp. This figure shows a lumped mass connected to 

the barrels via rigid elements. 

front of the thrust bearing. Instead of obtaining the continuovis d 3 mamic behavior 

j 

of the gun body, it has been broken up into two distinct parts with the thrust 
I bearing serving as the divider. The forward part of the gun body acts like a short 

dnim, anchored at the thrust bearing. By itself, its finite element representation is 
too coarse and leads to unrealistic modeUng behavior. This behavior appears in 
the higher order modes of vibration. If the forward gun body is coupled to the rear 
gun body, as it should be, then its finite element representation is not too coarse 
and no further refinements are needed. 
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How is this continuity achieved? The brick elements used to represent the 
thrust bearing should be replaced by beam elements. In I-DEAS, a beam element’s 
cross-section can be defined in as much detail as desired. The inner and outer 
race, as shown in Figure 3-15, could be defined as a beam cross-section. By 
replacing the bricks with beams, the degrees of fi'eedom are then matched with the 
thin shells, since beams also have six degrees of freedom. Therefore, the rear gtm 
body and forward gun body are now coupled and continuous dynamic behavior 
would be achieved. 

In making this transition to beam elements, the springs used to model 
the ball bearings would also have to be modified. Now rigid elements would be 
used to extend firom the centroid of the inner and outer races to a common, 
coincident point between them. The springs wotdd be co-located, one for radial 
stiffness, axial stiffness and torsional stiffiiess. This refinement would also serve 
the purpose of removing any moment arms that currently exist in the model due to 
the finite distance the springs extend between the inner and outer race. 

More attention could be used when modeling the details of the gun body 
housing. Currently, it is represented as a cylinder, whereas in reality, there are 
many intricate details that might be important to model more accurately. For 
example, the ammunition chuting and a hydraxilic motor coupled to the gun drive 
could be better studied now that the gun is on-sight. 

3. Rotor 

The rotor is a cylindrical component represented with brick elements as 
shown in Figure 3-16. The dynamics of the rotor are not critical to understanding 
the gun, but could be better represented by a lumped mass and inertia properties 
while keeping a portion of the brick elements for use in attaching to the thrust and 
needle bearings. In this way, the degrees of freedom of the model are reduced, 
while preserving the load path from the excitation force to ground. The bricks 
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would be important if there was some concern that shear forces were causing it to 
deform, or if thermodynamic properties were being modeled. Since neither of these 
characteristics are being modeled, the use of brick elements for the whole rotor is 
unnecessary. The portions used to couple the barrel assembly to the bearings 
should be kept, and should be connected via rigid elements to a lump)ed mass at its 
center of gravity. 

4. Stub Rotor 

The stub rotor, constructed much like the muzzle clamp as Figure 3-11 
shows, with webbing between two plates, could be further improved by dividing the 
single row of bricks for each plate into two rows. The stub rotor is a critical 
component in the dynamics of the gun, perhaps more so than any other component. 
Its rigidity is important in determining how much the gun barrels will displace and 
has been experimentally studied by Martin Marietta to calculate its characteristics 
[Ref. 16]. The component should be modeled using more elements instead of fewer 
elements. By including two rows of elements, its bending behavior can be better 
represented. One aspect of the stub rotor and barrel clamps that is diffi cult to 
model is the nonlinear effects that gaps play in the gun dynamics. The current 
version of I-DEAS does not allow this option, but a new release will allow the 
nonlinear gaps to be modeled. 

5. Recoil Adapters 

The recoil adapters have been the most challenging components to model. 
The current model defined them in terms of "ears” that stood out on either side of 
the gun body that were connected to ground via rigid elements and recoil springs 
as shown in Figure 3-16. The "ears" were solely used to represent the adapter’s 
mass and to model the correct geometric position of the attachments. The model 
could be simplified if the "ears" were replaced with lumped mass elements. In 
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addition, the rigid elements could connect the gun body to the springs directly 
while maintaining the correct offset of the pivot points. These changes are needed 
in order to provide better coupling of the recoil springs to the gun body. As 
discussed earlier, problems occur when there is a discontinuity between degrees of 
freedom of elements. In this case, it occurs between the thin shells of the body to 
the bricks of the "ears." This recommendation is better than trying to accurately 
model the recoil adapters. It is likely that little would be gained for such detailed 
modeling and would significantly add to the model’s complexity. 

B. VALIDATION OF THE MODEL THROUGH TESTING 

After constructing the finite element model of the gun and obtaining its 
normal modes, the next logical step would be to compare both the mode shapes and 
their associated fipequency with those of the real gun. However, it was only until 
late in this study that a M61A1 PHALANX gun was received at the school. After 
proper mounting, a modal analysis of the gun could then be performed and the 
results compared to the model. The equipment already exists on-sight for this type 
of testing. In addition to testing the whole gtm as one unit, it could be broken 
down into a barrel cluster and gun body which would then be tested separately. In 
this way, discrepancies could be localized to determine which components needed 
to be better modeled. 

Another comparison that can be made using the real gun is that of finding the 
proper damping values. Although the damping value of the recoil adapters was 
provided by Martin Marietta and verified through analytical means in Chapter VI, 
it would be useful to obtain these values experimentally. It is uncertain whether 
the damping values vary much from one adapter to another, or the effect wear 
plays on them. More importantly though, all other normal modes of the gun were 
estimated to have a damping value of two percent. The damping values of the gun 
model could be correlated with the real gun during its modal analysis. 
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As mentioned earlier, it is also uncertain how much free play exists in the 
gun and the amount each component contributes to the total. This information 
could also be determined experimentally. Individual components of the gun should 
be weighed in order to provide another means of comparison so that the modeling 
could be improved. 

C. ADDITIONAL MODELING 

1. Modifying the Existing Firing Force 

In an efifort to get through a complete modeling cycle, there were many 
variations of the existing firing force that were put off until a later time. One of 
the modifications which would certainly add to the accxnracy of the model is to 
ramp up the firing force. The real gun does not immediately jump to a firing rate 
of 4500 shots per minute. This could be accomplished by spacing out the firing 
impulses for the first ten to twenty shots according to data available from the 
Naval Surface Warfare Center in Dahlgren or Port Hueneme. Further 
modifications would be to model the difierent firing rates themselves. This study 
only looked at the firing rate of 4500 shots per minute. However, the existing 
PHALANX gun can select a rate of fire of either 3000 or 4500 shots per minute. 
There has been a proposal to investigate a firing rate of 6000 shots per minute. If 
the mass of the penetrator were to change or if the penetrator’s velocity were 
modified, the firing force would change and could easily be mod''led. The gun is 
capable of longer engagements and longer bursts of shots. Longer firing bursts 
shoiild be modeled. An additional and perhaps important firing force to model 
would be that of an occasional misfire. In this case, the g\in has a chance to 
translate forward after the misfire and then be slammed back as the next round is 
fired. 
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2. Adding Additional Forces 

In this study, only the firing force was represented. The impulse of 
loading each round may also affect the gun. In discussions with Martin Marietta, 
this force is not insignificant since the round has a tight fit in the barrel. This 
loading force is applied as the barrel passes through its lowest point of circular 
motion. While the fi:equency of the loading force must be the same as the firing 
rate, the phase relation may allow it to add to or subtract from the firing force. 

The effects of rotation can also be modeled in I-DEAS. Here the damping 
matrix is modified to represent the effect of rotation. Even though the modeled 
gun is not rotating, it can be "tricked” into producing the same dynamic results as 
if it were rotating by modifying its damping matrix. This involves adding off- 
diagonal terms in the viscous damping matrix. The normal modes of vibration are 
then resolved as before, and the excitation forces applied. This technique is only 
good for one firequency, and the model will have to be resolved for each different 
firing rate. 

Another force to model is the variation in rotation of the gim by adding 
an impulse perpendicular to its axial length. Even though this force may be small, 
its effect could be large since it might be a contributing perturbation to significant 
dynamic forces such as whirl. 

3. Additional Modeling 

With the current PHALANX gun model partially validated, it is possible 
to move on and use it as a tool. The U. S. Navy is interested in the possibility of 
reducing the PHALANX gun’s dispersion by constructing a gun barrel restraint like 
the one shown in Figure 7-2. This structure wiU be clamped at its base to the 
isogrid of the gun mouint. The other end will use a radial bearing to clamp onto 
the gun barrels. Its purpose is to restrict the most extreme motions of the barrel 
cluster. A finite element analysis has not been used to date to study the motion of 


123 





thin structure and optimize its design. The blueprints have been received, and the 
Naval Postgraduate School is poised to take on this task [Ref. 24]. The computer 
model will be able to study the mode shapes of this combined system and check if 
there is any problem resonant frequencies. By studying the mode shapes, it wUl be 
possible to determine the best location to clamp the bearings. This system can be 
driven and damped to calculate the new barrel tip displacements. The muzzle 
restraint structure could be modified in various ways to find its optimum design 
and placement. 



for Strength 


Figure 7-2. Proposed assembly of the PHALANX muzzle restraint [Ref. 24 

Another modeling task is to swap out the cxirrent set of finite element 
barrels in exchange for the proposed optimized barrels. This will help determine 
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whether the new barrels will pitch more than the old ones. The new barrels will be 
much heavier, nearly 18 kg (39 lb) compared to the current 8 kg (18 lb) barrels. 
They will be much longer, roughly 200 cm (79 in) compared to the cxirrent 152 cm 
(59 in) barrel length. The effect of this extra weight and length on the gun mount 
and its dispersion has not been modeled to date. The NPS model would be an 
excellent tool for making this analysis. The firing force would have to be modified 
to take into account the extra recoil force associated with the heavier enhanced- 
lethality projectile designed to be used by these barrels. It has a greater mass of 
0.290 kg (0.64 lb) than the old projectile mass of 0.245 kg (0.54 kg), and will be 
fired at a greater velocity increasing the specific impulse on the gun mount. An 
additional modeling step would be to use the proposed muzzle restraint, as 
described above, with these new barrels and solve for its dynamic response. 
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VIII. CONCLUSIONS 


Using a finite element model to understand the complex dynamics of the 
PHALANX gun has made possible new areas of analysis and component design. 
This has been the first effort to use the sophisticated I-DEAS computer code to 
study the gun. According to the Naval Surface Warfare Centers in Dahlgren and 
Port Hueneme, it has also been the most "technologically advanced" modeling 
research of this important weapon system [Ref. 25]. Not only does it provide a 
thorough and detailed model of the gun where various driving forces can be applied 
and the barrel tip response calcxilated, but it is a working tool to assist the 
PHALANX program. In particvilar, each individual component of the gun can be 
studied in detail to understand how much it affects the overall dispersion. This is 
an especially important tool today when budget hmitations allow only modest 
improvements to be made to the gun. Furthermore, new components can be tested 
and their performance predicted before any actual hardware tests. This can be 
applied to the proposed muzzle restraint and set of optimized barrels, before any 
prototypes are ever constructed. The cost saving and, not to mention, time saving 
aspects of this method provides an invaluable resource to the program. This t}^ 
of research may not have been carried out before by private industry or the U. S. 
Navy because the computer code is expensive and has only been made easily 
accessible in the past few years. The Naval Postgraduate School receives an 
educational discount to obtain an on-sight hcense for the code. In addition, the 
high concentration of professors with doctorates and U. S. Navy officer/students 
makes fur a unique and productive environment. The school has been proposed as 
the center of all future modeling research of the PHALANX gun [Ref 26]. In 
summary, the U. S. Navy is now using state-of-the-art technology to understand 
and improve old technology in a cost saving way. 
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The primary goal of thia study has been to construct an accurate finite 
element model of the PHALANX gun and demonstrate the potential of such a 
sophisticated computer code. This goal has been achieved. Parts of the model 
have been validated with the real gun, and the results of this study have been 
presented to the several organizations involved in improving the PHALANX gun 
[Ref 21]. The model is poised to take on most of the new design changes that the 
group proposes. The model is now considered an integral part of the decision 
making process. With an actual g\in on-sight, the model can be completely verified. 

Although the primary purpose of constructing this finite element gun model 
has been achieved, some preliminary analysis was conducted along the way. One 
of the first areas of concern to be noticed was the gun’s normal mode of vib'-ation 
around 14 Hz. This first flexible mode with associated pitching motion would 
contribute to dispersion, if it is sufficiently excited. The gun rotates at a frequency 
of 12.5 Hz making this rotational force a contender for exciting this mode. 
Frequency response functions taken of the barrel tip indicates that this mode is 
being excited by the principal driving force, of the gun firing. This discovery was 
somewhat unexpected, but upon analysis of the moment arm created from the 
driving force and the gun moimting positions, it is reasonable that it is excited. As 
the giui fires, the moment arm causes the gun to rock back and pitch upward 
coupling the firing recoil into other planes of motion. Fourier transforms performed 
of actual dispersion tests appear to correlate this 14 Hz normal mode [Ref 22). 

As far as design analysis of the g\in, the most critical component contributing 
to dispersion was revealed by this study as the double angular contact bearing, or 
thrust bearing. When various bearing stiffiaess values were tested, the effect on 
barrel tip displacement was significant. Confirmation of this fact came from the 
field where the dispersion pattern improved tremendously as a worn thrust bearing 
was replaced with a new one [Ref. 23]. No data is available on the rate of wear of 
the thrust bearing. In fact, the bearing itself is not inspected while on ships, but 
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only when the gun is returned for depot level maintenance. The faci that the 
bearing is not pre-loaded is a leading cause of the physical ejQTect known as whirl 
that occurs on an overhung shall as in the PHALANX gun. Dispersion data also 
suggests the possibihty that whirl is occurring. The addition of the new, heavier 
and longer barrels may significantly aggravate whirl and counteract any dispersion 
improvements achieved by the new barrels, especially if the bearing is not 
redesigned to support these new barreb. Unfortunately, this appears to be the 
case. Not only would the added torque of the new barrels lead to a greater amount 
of whirl, but it would probably cause a greater amount of wear on the bearing 
which would further worsen the problem. 

One of the first validation tests of the gun model was to drive it with the 
principal driving force of the rounds firing. In the response plots jhown in Chapter 
VI, the firing gun barrel acted as predicted. That is, the gun set back the observed 
1.27 cm and oscillated about this new equilibrium point. Barrel tip displacement 
in the pitch and yaw directions corresponds well with the observed 2 milliradian 
dispersion. Elliptical dispersion patterns were also predicted for the first time. 

In this first finite element study of the PHALANX gun, the potential for its 
future use has been documented in this effort. It is a cost effective tool that has 
been designed to help the U. S. Navy not only understand the current gun’s 
complex dynamics, but those of any future g\in as well. As a ship’s last layer of 
defense, the PHALANX gun provides a valuable function for the ship and her crew. 
This is especially true in the near future where the U. S. Navy has extended its 
operating envelope into new environments, where the l^t layer of defense may be 
the first layer. 
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APPENDIX A. GUN BARREL NORMAL MODES 


The following graphs are experimental data of the normal modes of a single 
PHALANX gun barrel. The barrel was suspended by two low stiffness cords and 
tapped with an instrumented hammer. The resulting vibrational information was 
recorded with a transducer mounted on the opposite side of the hammer impact 
point. The normal mode frequencies were determined through 3000 Hz and appear 
as spikes in the following graphs. 
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APPENDIX B. GUN MODEL MODE SHAPES 


The following pictures are snapshots of the mode shapes used in the dynamic 
analysis of the gun. There are 20 flexible modes listed in ascending order by 
frequency. 
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